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PREDICTING  THE  PERFORMANCE  OF  BAND  7 
COMMUNICATION  SYSTEMS 
USING  ELECTRONIC  COMPUTERS 


Donald  L.  Lucas  and  George  W.  Haydon 


Radio  system  parameters  are  combined  with 
geophysical  and  ionospheric  characteristics  to  pre¬ 
dict  the  performance  of  high  frequency  sky-wave 
communication  circuits  through  the  use  of  electronic 
computers.  A  program  is  presented  to  compute 
Maximum  Usable  Frequencies,  Optimum  Traffic 
Frequencies,  Lowest  Useful  Frequencies ,  probable 
mode  of  propagation,  angle  of  arrival,  circuit 
reliability,  system  loss,  available  s ignal-to -noise 
and  field  strength.  Numerical  representation  is 
used  for  all  parameters  not  expressed  in  closed 
mathematical  form. 

I.  INTRODUCTION 

Long  distance  high  frequency  radio  systems  are  subject  to 
marked  variations  in  performance,  most  of  which  are  directly  related 
to  changes  in  the  ionosphere.  Changes  in  the  ionosphere  result  in 
variations  in  the  maximum  frequency  which  will  be  returned  back  to 
earth,  in  the  strength  of  the  radio  waves  due  to  increased  absorption, 
in  optimum  vertical  angles  of  wave  arrival  and  departure,  and  in  the 
background  noise  level  as  atmospheric  radio  noise  is  propagated  from 
great  distances.  Effective  operation  of  long  distance  high  frequency 
radio  systems  increases  in  proportion  to  the  ability  to  predict  variations 
in  circuit  performance  and  thereby  optimize  frequency  selection,  antenna 
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choice  and  other  circuit  parameters  to  capitalize  on  anticipated  ionos¬ 
pheric  conditions.  It  is  the  purpose  of  this  report  to  present  a  current 
computer  routine  which  has  been  developed  to  combine  the  more  pre¬ 
dictable  ionospheric  characteristics  with  circuit  parameters  to  calculate 
the  expected  performance  of  high  frequency  sky-wave  radio  systems. 
Basic  ionospheric  characteristics  which  are  considered  predictable 
enough  to  be  useful  are: 

(1)  The  monthly  median  of  the  ordinary  ray  vertical  incident 
critical  frequencies  of  the  F2  layer  (foF2). 

(2)  The  monthly  median  of  the  factors  relating  the  vertical 
incident  critical  frequency  to  Maximum  Usable  Frequency  on  a  3000  km 
path  (M-3000  factor). 

(3)  The  monthly  median  minimum  virtual  height  of  the  F2  layer. 

The  following  ionospheric  characteristics  are  derived  in  the 

prediction  process: 

(1)  The  E-layer  critical  frequency  as  a  function  of  solar  activity 
and  angle  of  the  sun. 

(2)  Ionospheric  absorption  as  a  function  of  operating  frequency, 
gyro  frequency,  solar  activity,  angle  of  the  sun,  and  vertical  angle  of 
departure . 

(3)  Factors  relating  the  E-layer  critical  frequencies  at  different 
departure  angles. 

(4)  Factors  relating  F-layer  critical  frequencie s  at  different 
departure  angles. 

The  following  geophysical  data  are  utilized: 

(1)  E-region  gyro  frequency. 

(2)  F -region  gyro  frequency. 

(3)  Geomagnetic  latitudes. 

(4)  World-wide  distribution  of  land  and  sea. 
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(5)  World-wide  distribution  of  atmospheric  noise  including 
seasonal,  diurnal  and  frequency  variations. 

(6)  Cosmic  noise  including  frequency  dependence. 

(7)  Man-made  noise  in  the  receiving  area  including  frequency 
dependence. 

The  following  circuit  data  are  needed: 

(1)  Transmitter  location. 

(2)  Receiver  location. 

(3)  Transmitter  power. 

(4)  Transmitting  antenna  (physical  dimensions  and  orientation). 

(5)  Receiving  antenna  (physical  dimensions  and  orientation). 

(6)  Type  of  traffic  (modulation,  type  of  intelligence,  speed  of 
transmission,  required  quality). 

The  above  data  are  combined  to  determine  the  following  perfor¬ 
mance  characteristics  of  the  circuit: 

(1)  The  Maximum  Usable  Frequency  (MUF). 

(2)  The  Optimum  Traffic  Frequency  (FOT). 

(3)  System  loss  (ratio  of  power  received  to  power  transmitted). 

(4)  Field  strength  at  the  receiver. 

(5)  Monthly  median  of  hourly  median  available  signal -to-noise 

ratio. 

(6)  Circuit  reliability  (per  cent  of  days  the  hourly  median  signal- 
to-noise  ratio  can  be  expected  to  equal  or  exceed  a  given  value). 

(7)  The  Lowest  Useful  Frequency  (LUF),  the  frequency  below 
which  the  circuit  reliability  is  expected  to  be  less  than  90%. 

In  addition  to  the  above  circuit  performance  characteristics,  the 
following  other  circuit  characteristics  are  determined: 

(1)  The  short  great  circle  distance. 

(2)  Bearing  of  the  transmitter  at  the  receiver. 
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(3)  Bearing  of  the  receiver  at  the  transmitter. 

(4)  Propagation  path  having  the  lowest  system  loss,  e.g.  ,  2  hops 
via  F  layer  =  2F. 

(5)  Vertical  reception  angle  associated  with  the  best  propagation 

path. 

The  solution  of  the  problem  is  divided  into  two  parts:  (1)  an 
estimation  of  the  available  signal,  and  (2)  an  estimation  of  the  required 
signal.  The  available  signal  depends  upon: 

(1)  The  transmitter  power. 

(2)  The  transmitter  antenna  gain. 

(3)  Loss  due  to  spreading  of  the  radio  energy  as  it  propagates 
to  greater  distance. 

(4)  Losses  in  the  ionosphere  due  to  absorption. 

(5)  Losses  at  each  ground  reflection. 

(6)  The  gain  of  the  receiving  antenna. 

The  required  signal  depends  upon: 

(1)  The  atmospheric  noise  at  the  receiver  location. 

(2)  The  man-made  noise  at  the  receiver  location. 

(3)  Cosmic  noise. 

(4)  The  required  signal -to -noise  ratio  depending  Upon  type  and 
quality  of  service  desired. 

The  most  basic  calculation  in  the  estimation  of  available  signal 
power  is  the  system  loss  calculation.  The  ionospheric,  geophysical 
and  circuit  characteristics  are  combined  to  estimate  a  quasi-minimum 
system  loss  (the  minimum  loss  expected  on  any  day  of  the  month  at  the 
hour  for  which  the  calculations  are  made).  The  distribution  of  losses 
above  the  quasi-minimum  have  been  empirically  determined  and  are  a 
function  of  geographic  location  and  time  of  day.  These  statistical 
distributions  are  used  to  account  for  the  numerous  other  factors  which 
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contribute  to  variations  between  predicted  signal  levels  and  those  ob¬ 
served  on  a  given  day.  These  factors  include  polarization  mismatch 
between  the  signal  and  the  antenna,  focusing  by  the  ionosphere,  defo- 
cusing  by  the  earth,  variations  between  theoretical  and  actual  antenna 
performance  and  day  to  day  variations  in  ionospheric  layer  heights, 
ionospheric  absorption  and  critical  frequencies. 

Basic  F2-layer  data  for  long  term  predictions  are  numerical 
maps  of  the  F2  region  as  determined  from  world-wide  observations 
using  1954  as  typical  of  a  low  solar  activity  period  and  19  58  as  typical 
of  a  high  solar  activity  period.  In  addition  to  19  54  and  19  58,  F2  data 
from  other  periods  of  high  and  low  solar  activity  have  been  incorporated 
in  the  March,  June,  September  and  December  predictions,  since  these 
seasonal  extremes  are  often  used  in  long  term  planning. 

Numerical  maps  of  the  F2  layer  for  monthly  predictions  will  be 
a  part  of  the  "Ionospheric  Predictions"  CRPL  Series  D  issued  by  the 
Central  Radio  Propagation  Laboratory  starting  with  the  January  1963 
issue. 

Sporadic  E  and  Fl-layer  predictions  are  considered  in  this  report, 
only  insofar  as  propagation  by  these  modes  enter  into  the  empirically 
determined  signal  distributions. 

II.  PROGRAM  DESCRIPTION 

The  computer  routine  is  based  on  established  manual  methods 
and  assumes  a  working  knowledge  of  these  methods  [  Laitinen,  1949], 
[Haydon,  1962],  [NBS  Circular  No.  462],  along  with  a  familiarity  of 
Fortran  II  computer  language  [McCracken,  1962], 

The  computer  method  described  is  designed  for  any  (32K  IBM- 
7090  class)  computer  with  the  necessary  compiler  and  tape  units.  The 
Fortran  II  computer  language  assures  flexibility  and  ease  of  modification 
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as  better  knowledge  of  the  various  parameters  becomes  available.  The 
program  provides  the  communicator  with  a  completely  general  method 
for  predicting  the  performance  of  any  sky-wave  high  frequency  radio 
system.  The  system  performance  calculations  may  be  performed  for 
any  month  of  the  year,  and  any  degree  of  solar  activity  for  any  hour  of 
the  day.  Numerical  mapping  is  included  for  those  parameters  not 
evaluated  directly  by  the  closed  mathematical  expressions  (Section  IX). 

The  program  is  designed  to  be  used  with  either  a  "long-term" 
prediction  data  tape  which  utilizes  coefficients  foF2  and  M-3000  factors 
representing  periods  for  high  and  low  solar  extremes,  or  a  "short-term 
data  tape  which  utilizes  monthly  coefficients  issued  by  CRPL  approx¬ 
imately  three  months  in  advance. 

1.  MUF-FQT 

The  Maximum  Usable  Frequency  and  Optimum  Traffic  Frequency 
are  based  on  F2-region  numerical  mapping  of  ionospheric  character¬ 
istics  recently  developed  by  CRPL  [Jones  and  Gallet,  1961],  and  a 
semi-empirical  relationship  between  the  sun' s  zenith  angle  and  solar 
activity  for  the  regular  E  region  [Haydon,  1962].  Conventional  pre¬ 
diction  methods  [NBS  Circular  No.  462]  using  control  points  along  the 
great  circle  2000  kilometers  from  each  terminal  for  the  F2  layer,  and 
control  points  1000  kilometers  from  each  terminal  for  the  regular  E 
layer  are  used  for  paths  greater  than  4000  kilometers.  The  midpoint 
of  the  path  is  considered  as  the  control  point  for  both  E  and  F  layer  on 
paths  less  than  4000  km. 

Since  the  basic  data  for  the  predictions  are  monthly  median  values 
predictions  for  specific  days  within  the  month  are  not  available. 
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2.  System  Loss 

1.  System  loss  computations  are  based  on  ray-path,  theory  with 
path  limitations  based  on  ray-path  geometry  between  a  spherical  earth 
and  a  concentric  ionosphere.  Limitations  fall  into  two  categories, 
arrival  angle  and  penetration  frequency,  and  the  probable  paths  for  a 
given  circuit  at  a  given  hour  and  frequency  are  chosen  on  the  basis  of 
these  limitations. 

F2-region  propagation  is  considered  unlikely  in  these  computations 
if  the  angle  at  the  F2  region  is  too  sharp  for  F2-layer  support  or  if  the 
penetration  frequency  of  the  E  region  rises  to  a  value  which  will  support 
the  operating  frequency.  In  this  last  case,  E  region  transmission  is 
geometrically  possible.  Multiple  hop  transmission  for  a  given  circuit 
is  possible  solely  by  the  E  region  or  the  F2  region,  but  the  computer 
routine  also  includes  the  case  where  some  hops  of  a  given  path  are 
supported  by  the  E  region  and  other  hops  supported  by  the  F2  region. 

Three  F2-region  paths,  two  E-region  paths,  and  one  E-F2-region 
path  are  inspected  for  a  given  circuit  at  a  given  hour  for  a  given  fre¬ 
quency.  The  path  that  is  geometrically  possible  with  the  least  theo¬ 
retical  loss  is  assumed  to  be  the  most  useful  mode  with  a  loss  typical 
for  the  circuit.  The  total  field  is  not  included  in  this  routine. 

2.  Layer  heights  used  for  the  F2  region  are  minimum  vertical 
heights  averaged  for  a  given  month.  E-layer  height  is  assumed  to  be 
110  km.  The  average  height  of  the  absorbing  region  is  assumed  to  be 
100  km. 

3.  Ionospheric  layer  heights  at  the  reflection  points  are  averaged 
in  calculation  of  the  vertical  radiation  angle  (A)  and  the  angle  of  incidence 
at  the  ionosphere  (<f>').  Unequal  take-off  angles  in  multi-hop  cases  are 
not  included. 
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4.  Absorption  in  the  lower  regions  is  represented  by  a  semi- 
empirical  relationship  involving  the  angle  of  incidence  at  the  ionosphere, 
the  gyro  frequency  in  the  absorbing  region,  the  zenith  angle  of  the  sun, 
solar  activity,  and  the  operating  frequency. 

5.  Estimation  of  the  ground  reflection  losses  is  evaluated  for 
either  land  or  sea  reflections.  Shore  line  and  mixed  paths  are  not 
included  and  random  polarization  of  the  downcoming  waves  is  assumed. 

6.  Convergence  and  divergence  of  the  ionospheric  waves  are 
included  in  the  routine  only  in  the  empirical  determined  system  loss 
distributions . 

7.  Theoretical  vertical  and  horizontal  reflection  coefficients 
are  evaluated  for  the  ground  losses  and  the  antenna  patterns. 

8.  Receiving  antenna  response  (gain)  is  assumed  to  be  equal  to 
the  antenna  gain  of  an  identical  antenna  used  as  a  transmitting  antenna. 
All  gains  are  relative  to  an  isotropic  in  free  space. 

The  complete  computer  routine  contains  a  subroutine  for  three 
general  types  of  antennas. 

a.  Terminated  rhombics  at  any  height  above  the  terrain,  any 
leg  length  and  tilt  angle.  Poor  earth  (O’  =  .  001  mhos/  meter;  €  =  4)  is 
assumed  as  the  terrain,  but  any  ground  conductivity  and  dielectric 
constant  could  be  used.  Off-azimuth  radiation  may  be  computed  for 
the  rhombics. 

b.  Horizontal  half-wave  dipoles  for  any  height  above  finite  earth. 

c.  Vertical  ship  antennas  for  constant  length  elements  or  for 
multiples  of  any  wave  lengths  not  less  than  .  02  Practical  efficiencies 
of  grounded  verticals  less  than  one -quarter  wave  length  are  computed 
and  included  in  the  antenna  gain. 
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9.  Signal  distribution  is  assumed  to  vary  with  geomagnetic 
latitude  and  length  of  the  circuit  and  quasi-minimum  losses  are 
adjusted  to  median  losses  in  terms  of  these  empirical  distributions 
[Haydon,  1962], 

10.  Off-great-circle  transmissions  are  not  included  in  these 
calculations . 

3 .  Field  Strength  Option 

The  field  strength  option  is  computed  similar  to  the  system  loss, 
then  mathematically  related  to  the  system  loss  (Section  IX). 

4.  Radio  Noise 

The  CCIR  world-wide  maps  of  atmospheric  noise  are  used  as 
representative  of  the  distribution  of  noise  [CCIR  Report  No.  65], 

Man-made  noise  estimates  are  based  upon  the  type  of  receiving 
area  [Haydon,  1962].  Measured  values  of  man-made  noise  are  also 
acceptable. 

Galactic  noise  at  the  receiver  site  [  CCIR  Report  No.  65]  is 
estimated  by  a  least  squares  polynomial  in  frequency.  The  highest  of 
the  three  noise  types  is  assumed  to  predominate  and  the  total  noise 
field  is  not  estimated.  All  noise  values  are  expressed  in  db  relative 
to  1  watt  for  a  one -cycle  bandwidth.  The  1  Me  /  s  atmospheric  noise 
above  KT  B  is  adjusted  to  a  1  c  /  s  band  relative  to  one  watt  by  least 
squares  polynomial  representation  of  the  seasonal  maps  [Lucas  and 
Harper,  1962], 

5.  Available  Signal-To -Noise 

The  monthly  median  of  the  hourly  median  available  signal  is 
compared  with  the  average  noise  in  a  four -hour  time  block  at  the 
receiver  site.  This  gives  an  estimate  of  the  monthly  median  of  hourly 
median  total  signal  pov/er  relative  to  noise  in  a  one -cycle  band. 
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6.  Required  Signal-To-Noise 

The  required  signal -to-noise  depends  upon  the  service  required, 
e.g.  ,  error  rate  in  teletype,  intelligibility  of  voice,  etc.  A  basic 
signal -to-noise  requirement  based  upon  a  steady  static  signal  and 
random  noise  adjusted  for  signal  fading  degradation,  diversity  improve¬ 
ment,  etc.  ,  determines  the  required  signal-to-noise. 

7.  Reliability 

Reliability  is  estimated  at  discrete  frequency  bands  by  calculating 
the  ratio  of  the  available  monthly  median  signal-to-noise  of  Paragraph 
5  to  the  required  hourly  median  signal-to-noise  of  Paragraph  6.  Using 
the  above  parameter  and  the  operating  frequency,  theoretical  reliability 
is  obtained  for  a  given  circuit  of  a  given  length  and  geomagnetic  latitude 
by  use  of  statistical  two-dimensional  maps. 

8.  Lowest  Useful  Frequency 

The  Lowest  Useful  Frequency  is  calculated  by  stepping  through 
the  HF  band  until  a  reliability  of  90%  or  more  is  obtained,  then  stepping 
down  in  small  increments  of  frequency  until  a  frequency  is  found  that 
has  (within  the  tolerance  of  the  frequency  increment)  90%  reliability. 
This  frequency  is  the  classical  90%  LUF. 

III.  PROGRAM  USE 

The  IBM- 7090  Fortran  II  listing  is  shown  in  Section  VIII.  The 
cards  corresponding  to  these  listings  with  the  associated  control  cards 
necessary  for  the  specific  monitor  being  used  are  all  that  is  necessary 
to  compile  the  program.  A  pictorial  view  of  a  sample  deck  ready  to 
run  under  the  IBM-7090  monitor  is  shown  in  Figure  17.  System  sub¬ 
routines  used  in  the  program  are  as  follows: 


Subroutine 
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Call 

Argument 

1. 

Sine 

SINF  (x) 

(radians ) 

2. 

Arc  sine 

ASINF  (x) 

(.  1  <  x  <  1) 

3. 

Cosine 

COSF  (x) 

(radians ) 

4. 

Arc  cosine 

ACOSF  (x) 

(-  1  <  x  <  1  ) 

5. 

Absolute  value 

ABSF  (x) 

(floating ) 

6. 

Absolute  value 

XABSF  (x) 

(fixed) 

7. 

Arc  tangent 

ATANF  (x) 

(  -  00  <  x  —  oo 

8. 

Square  root 

SQRTF  (x) 

8 

V! 

X 

VI 

o 

9. 

Transfer  of  sign 

SIGNF  (x,  y ) 

(arg  1,  arg  2) 

10. 

Logarithm  (Base  10) 

LOG  1  OF  (x) 

(0  <  x  <  oc  ) 

11. 

Natural  logarithm 

LOGF  (x) 

(0  <  X  2?  oc  ) 

12. 

Integer  function 

INTF  (x) 

(  _  00  <  X  ^  oc 

The  entire  program  including  common  variables  and  system 
subroutines  occupies  approximately  24000  decimal  locations. 

Routines  included  in  the  program  with  a  brief  description  of  their 
function  are: 

MUFLUF  --  (2550  decimal  locations) 

This  is  the  main  program.  It  calculates  the  Maximum  Usable 
Frequencies,  Optimum  Traffic  Frequencies,  and  other  parameters 
needed  not  only  in  the  calculation  of  these  critical  frequencies  but  also 
in  one  or  more  subroutines. 

LUFFY  --  (3857  decimal  locations) 

This  subroutine  is  larger  and  more  complex  than  MUFLUF,  but 
is  a  subprogram  because  the  option  is  asked  for  in  MUFLUF.  This 
routine  calculates  the  systems  performance  of  a  given  circuit,  i.  e.  , 
system  loss,  reliability,  signal-to-noise  and  field  strength.  The  sub¬ 
programs  that  follow  are  used  by  LUFFY  to  perform  a  specific  task 
needed  in  the  system  performance  calculations. 


-  12  - 


POLY  --  (107  decimal  locations) 

This  subroutine  is  called  upon  to  evaluate  all  parameters  that 
are  represented  by  Nth  degree  polynomials  of  the  power  series  in  X. 
Example: 


Y  =  (  A  x  °  +  A  x  1  +  •  .  .  +  A  x  n  Z 
0,0  0,1  0,  n  J 


i  0  1  A  1 

+  A  x  +  A  x  +  .  .  .  +  A  x  i  Z 
1,0  1,1  1,  n  J 


+  •  .  • 


+ 


+  A  1  x 
m,  1 


+  .  .  • 


+  A  x 
m,  n 


n  'N 
J 


Z 


m 


LOSS  --  (214  decimal  locations) 

This  subroutine  evaluates  the  ground  reflection  losses  as  a 
function  of  the  reflection  angle  and  operating  frequency  over  finite 
earth.  Random  polarization  is  assumed  in  the  calculations  as  the 
vertical  and  horizontal  reflection  coefficients  are  assumed  to  contri¬ 
bute  equally  in  the  polarization  effects  of  the  loss. 

This  routine  is  unique  in  that  complex  arithmetic  is  absent  in 
the  calculations  of  these  reflection  losses  [  Phillips,  1961], 

GAIN  --  (684  decimal  locations) 

4E 

This  routine  calculates  theoretical  power  gains  of  horizontal 
rhombics  in  three  dimensions,  horizontal  dipoles,  and  practical 
grounded  vertical  radiators,  all  over  finite  earth.  Practical  lower 
limits  of  gain  have  been  set  as  -10  db  relative  to  an  isotropic  in  free 
space. 

Cl  --  (411  decimal  locations) 

This  function  subroutine  evaluates  the  sine  and  cosine  integral 
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x 


s  in  v  dv 


v 


x 


Ci  (x)  =  0  I 


cos  v  dv 


v 


when  called  for  by  the  loss  and  gain  routines. 

NOISY  --  (215  decimal  locations) 

This  routine  calculates  the  world-wide  distribution  of  atmospheric 
noise  and  major  land  bodies  which  are  represented  by  Fourier 
coefficients. 

The  evaluation  is  made  at  a  given  longitude  and  latitude  of  the 
receiving  location  or  ground  reflection  point. 

CURVY  --  (1462  decimal  locations) 

This  routine,  although  not  pertinent  to  the  solution  of  any  problem 
encountered  in  system  performances,  produces  a  semi-line  graph 
representation  of  the  diurnal  variation  of  either  or  both  the  Optimum 
Traffic  Frequency  and  the  Lowest  Useful  Frequency. 

VREFCO  --  (230  decimal  locations) 

This  routine  calculates  the  ground  reflection  coefficient  (K^  )  for 
vertically  polarized  waves. 

NOTE: 

The  ground  reflection  coefficient  (K  )  for  horizontally  polarized 

ri 

waves  is  calculated  within  subroutine  GAIN. 


IV.  DESCRIPTION  OF  DATA  TAPES 


1.  "Long-term"  Data  Tape. 

The  binary  data  tape  for  the  storage  of  coefficients  representing 
those  variables  not  expressed  in  closed  form  has  a  packing  density  of 
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556  frames  per  inch.  Twelve  logical  records  -  one  for  each  month  - 
are  used.  Each  logical  record  is  composed  of  physical  records  256 
words  in  length.  Approximately  800  feet  of  tape  are  used. 

Each  logical  record  contains  six  one -dimensional' arrays ,  three 
three-dimensional  arrays,  one  two-dimensional  arrya  in  the  following 
order: 

IL  (4),  JL  (4),  KL  (4),  LK  (4),  JAL  (4),  Q  (20,60,4),- 

A  (10,  7,  14),  P  (29,16,6),  ABP(2,6) 

The  variables  IL,  JL,  KL,  LK,  and  JAL  are  limits  to  which  the 
Fourier  generation  of  the  foF2  and  M-3000  factors  are  parried,  e.g.  , 
IL  (1)  is  associated  with  the  two  dimensional  set  of  coefficients 

Q  (M.N,  1). 

The  "Q"  arrays  [Lucas,  1961]  contains  the  coefficients  used  in 
the  generation  of  the  foF2  and  M-3000  factors  in  the  following  order 
for  each  month: 


1.  Q  (M,N,  1)  foF2  -  low  solar  activity 

2.  Q  (M,N,2)  foF2  -  high  solar  activity 

3.  Q  (M,N,3)  M-3000  -  low  solar  activity 

4.  Q  (M,N,4)  M-3000  -  high  solar  activity 

The  "A"  array  contains  the  coefficients  needed  for  other  para¬ 
meters  not  represented  in  closed  form: 

A(M,N,1)  gyro-frequency  distribution  [  Lucas  ,  1961] 
A(M,N,2)  H'f2  layer  height  charts  [  Haydon,  1962] 

A  (M,N,3)  blank  to  leave  space  for  numerical  representation 
of  measured  patterns  of  antennas 
A  (M,N,4)  nighttime  frequency  dependence  of  atmospheric 
noise  [  Lucas  and  Harper,  1962]. 

A  (M,N,  5)  daytime  frequency  dependence  of  atrpospheric  noise 
A  (M,  N,  6)  man-made  noise  level  in  industrial  area  [Haydon, 

1  962] 
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A  (M,N,7)  man-made  noise  level  in  residential  area 
A  ( M ,  N ,  8 )  man-made  noise  level  in  rural  area 
A  (M,  N,  9)  man-made  noise  level  in  remote  unpopulous  area 
A  (M,  N,  10)  daytime  distribution  of  circuit  reliability  in 
temperate  regions  [  Haydon,  1962] 

A  (M,N,  11)  nighttime  distribution  of  circuit  reliability  in 
temperate  regions 

A  (M,  N,  12)  distribution  of  circuit  reliability  in  polar  region 
A  (M,N,  13)  distribution  of  short  circuit  reliability  in 
auroral  region 

A  (M,N,  14)  distribution  of  long  circuit  reliability  in 
auroral  region 

The  "P"  array  of  each  logical  record  contains  five  maps  of  the 
world-wide  distribution  of  atmospheric  noise  [Lucas  and  Harper,  1962] 
for  the  month  being  considered.  The  major  land  bodies  of  the  world 
are  also  represented  in  this  array. 


P  (M,N,  1) 

2000-0400  LMT 

P  (M,N,2) 

1600-2000  LMT 

P  (M,N,3) 

1200-1600  LMT 

P  (M,  N,  4) 

0800-1200  LMT 

P  (M,  N,  5) 

0400-0800  LMT 

P  (M.N.6) 

major  land  bodies  of  the  world 

ABP"  contains 

the  convergence  factors  for  the 

generation  of  the  "P"  matrix  [  Lucas  and  Harper,  1962]. 

The  elements  of  the  matrices  that  are  not  sensitive  to  month  are 


identical  for  all  logical  records. 

2.  "Short-term"  Data  Tape. 

(a)  Description  of  "Short-term"  Data  Tape 
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The  "short-term''  data  tape  is  identical  to  the  "long-term"  data 
tape  with  the  exception  of  the  £oF2  and  the  M-3000  coefficients.  On 
the  "short-term"  tape  each  particular  month's  coefficients  are  stored 
in  positions  used  by  the  coefficients  representing  high  and  low  solar 
activity  extremes  on  the  "long  term"  tape.  The  solar  activity  index 
associated  with  the  "short-term"  coefficients  must  be  used  when  pre¬ 
dictions  are  made  using  these  coefficients.  The  routine  which  follows 
is  designed  to  generate  and  update  a  "short-term"  data  tape  as  the 
coefficients  become  available.  It  is  designed  to  use  the  "long-term" 
or  "short-term"  tape  (Logical  #2  and  Logical  tape  =3) which  will  become 
the  new  "short-term"  data  tape.  Precede  the  data  deck  with  one  card 
which  contains  the  first  and  last  months  for  which  the  "short-term" 
tape  is  to  be  updated.  A  minus  1  card  punched  in  columns  55-56  must 
follow  each  individual  foF2  or  M-3000  coefficient  deck.  A  sample  deck 
for  one  month’ s  data  is: 


6£J0j6  ‘  4503  139  10  6  43 


37.0 


3. 0433946E  00 


L  t u  ff — T’rr  i"  rmirai  m  i 

|  ‘  63TX  4500  139  30  S  56  10  24  1  .37.0  0  0  ?.98639{55E  00 

I  6  6  M 


Cx; 


FORTRAN  STATEMENT 
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(b)  Program  to  Use  CRPL  "Ionospheric  Predictions"  Coefficients 
to  Generate  and  Update  "Short-term"  tape. 

C  ALWAYS  USE  MOST  RECENT  TAPE  FOR  UPDATING  SHORT-TERM  TAPE 

DIMENSION  0(20*60*4) ,  IL (4) *JL( 4) *KL(4) *LK ( 4 ) * JAL ( 4 ) *A(10*7*14)»P 
129*16*6) *ABP (2*6 ) 

1  FORMAT  (19X*2I4*4X*2I4*12X*2I5*E17*7) 

2  FORMAT  (51X»2I5*E17.7) 

101  FORMAT (212) 

READ  INPUT  TAPE  5*101  * MONTH 1 *MONTH2 

REWIND  2 

REWIND  3 

MONT H3= MONTH 2+1 

LOCK  =  MONTH 1-1 

I F ( LOCK )  4*55*4 

4  DO  3  11=1*  LOCK 

READ  TAPE  2*  IL,JL*KL,LK*JAL*Q*A*P1*ABP 

3  WRITETAPE  3 »  IL * JL *KL * LK . JAL *Q *  A *P 1  * ABP 
55  DO  100  1 1 = MONTH 1  * MONTH2 

READ  TAPE  2*  I L * JL *KL , LK, JAL *0* A *P 1  * ABP 
DO  600  KP=1 ,4800 
600  0 ( KP ) =0 # 

10=1 

5  READ  INPUT  TAPE  5*1*1 *J*K*L*M,N*P 
IL ( 10) =1+1 

JL ( 10) = J+l 
KL ( 1 0) =K+1 
LK ( 10) =L+1 
JAL ( 10 ) =2*J+1 
20  M=M+1 
N=N+1 

Q(M*N* 10) =P 

READ  INPUT  TAPE  5,2,M,N»P 
I F ( M )  40,20,20 
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40  10=10+2 

I F ( 4-10 )  6*5*5 

6  DO  7  10=2*4*2 

I L ( 1 0 )  =  I L ( IO-l ) 

JL(IO)=JL( IO-l) 

KL ( 10) =KL ( IO-l ) 

LK(IO)=LK( IO-l) 

JAL ( 10 ) = JAL ( IO-l ) 

DO  7  M  =  1 * 20 

DO  7  N=1 *60 

7  Q(M*N*IO)=Q(M*N*I 0-1 ) 

100  WRITETAPE  3*  IL* JL*KL*LK* JAL*0* A*P1 *ABP 
I F ( MONTH3-12  )  301  *301*401 
301  DO  201  K=M0NTH3 *  1 2 

READ  TAPE  2 * IL * JL *KL *LK * JAL *Q * A *P 1  * ABP 
201  WRITETAPE  3 ♦ IL * JL *KL *LK * JAL *0 * A *P 1 ♦ ABP 
401  REWIND  2 
REWIND  3 
CALL  EXIT 
END 

Minor  changes  in  the  Fortran  program  found  above  would  allow 
the  updating  of  the  "long-term"  tape  as  better  coefficients  become 
available . 

Extreme  caution  should  be  exercised  when  changes  are  being 
made  to  the  data  tape.  Fortran  statements  identical  to  the  read  state¬ 
ments  in  the  main  program  "MUFLUF"  should  be  used  to  re-write  the 
tape . 

Inspection  and  checking  of  coefficients  contained  on  the  tape  may¬ 
be  accomplished  by  merely  dumping  the  tape  off-line.  It  is  suggested, 
however,  that  for  easier  inspection  it  be  read  by  the  "read  tape"  state¬ 
ments  contained  in  "MUFLUF". 
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Data  tapes  for  the  IBM-7090  containing  all  the  logical  records 
may  be  obtained  from  the  Radio  Systems  Division,  National  Bureau  of 
Standards,  Boulder,  Colorado. 

V.  PROGRAM  OPTIONS 

Options  included  in  the  program  aid  in  the  solution  of  a  wide  range 
of  systems  problems.  Options  designed  to  be  printed  on  standard 
9  by  1 1  paper  are: 

1.  MUF  and  FOT  (Figure  1) 

The  Maximum  Usable  Frequency  (MUF)  and  Optimum  Traffic 
Frequency  (FOT)  are  printed  three  circuit  months  per  page. 

2.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE,  AND  RELIABILITY 
(Figure  2) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  circuit 
reliability  are  printed  one  circuit  month  per  page.  MUF  -FOT1  s  are 
calculated  for  each  hour  of  GMT.  Mode,  angle  and  reliability  are 
calculated  only  for  even  hours  of  GMT. 

3.  FOT  and  LUF  (Figures  3  and  4) 

The  Optimum  Traffic  Frequency  and  the  Lowest  Useful  High 
Frequency  are  printed  three  circuit  months  per  page.  The  Optimum 
Traffic  Frequency  is  always  given  each  hour  of  GMT  while  the  commun¬ 
icator  has  a  choice  of  each  hour  or  every  even  hour  of  GMT  for  the 
Lowest  Useful  High  Frequency. 

4.  MUF,  FOT,  MODE,  ARR IVAL  ANGLE  AND  SYSTEM  LOSS 
(Figure  5) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  monthly 
median  of  hourly  median  system  loss  are  printed  one  circuit  month  per 
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page.  MUF  and  FOT1  s  are  calculated  for  all  hours  of  GMT  while  mode, 
angle,  and  system  loss  are  calculated  for  even  hours  of  GMT. 

5.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE  AND  FIELD  STRENGTH 
(Figure  6) 

The  Maximum  Usable  Frequency  and  Optimum  Traffic  Frequency, 
theoretical  angle  of  arrival,  probable  mode  of  propagation  and  monthly 
median  of  the  hourly  median  field  strength  are  printed  one  circuit  month 
per  page.  MUF  and  FOT' s  are  calculated  for  all  hours  of  GMT  while 
mode,  angle,  and  field  strength  are  calculated  for  even  hours  of  GMT. 

6.  MUF,  FOT,  MODE,  ARRIVAL  ANGLE  AND  SIGNAL-TO- 
NOISE  (Figure  7) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency, 
probable  mode  of  propagation,  theoretical  angle  of  arrival  and  monthly 
median  of  the  hourly  signal-to -noise  at  the  receiving  antenna  terminals 
are  printed  one  circuit  month  per  page.  MUF  and  FOT' s  are  calculated 
for  all  hours  of  GMT  while  mode,  angle,  and  signal-to-noise  are  cal¬ 
culated  for  even  hours  of  GMT. 

7.  MUF,  FOT  AND  RELIABILITY  (Figure  8) 

The  Maximum  Usable  Frequency,  Optimum  Traffic  Frequency 
and  circuit  reliability  are  calculated  every  even  hour  of1  GMT  and 
printed  two  circuit  months  per  page. 

8.  FOT  GRAPHS  (Figure  9) 

The  Optimum  Traffic  Frequencies  are  plotted  frequencie s  vs. 

GMT  time  on  a  nonlinear  frequency  scale  by  the  on-line  printer. 

9.  FOT  AND  LUF  GRAPHS  (Figure  10) 

The  Optimum  Traffic  Frequencies  and  Lowest  Useful  High 
Frequencies  are  plotted  frequencies  vs.  GMT  time  on  a  nonlinear 
frequency  scale  by  the  on-line  printer. 
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The  above  options  may  utilize  one  of  two  input  data  tapes.  Long¬ 
term  predictions  are  available  from  a  data  tape  with  numerical  coeffic¬ 
ients  for  the  months  representing  both  high  and  low  solar  activity. 

These  extremes  are  used  to  predict  MUF  associated  with  the 
intermediate  sunspot  number. 

Short-term  predictions  are  available  from  a  data  tape  with 
numerical  coefficients  for  approximately  three  months  in  advance.  No 
interpolation  is  made  when  using  this  data  tape. 


VI.  COMPUTER  INPUTS 
1.  Circuit  Card  Format 


A  data  sheet  should  be  prepared  for  each  circuit  before  punching 
the  card.  The  number  in  the  left  hand  column  below  refers  to  infor¬ 
mation  on  the  data  sheet  (Figure  11).  The  fields  should  be  punched  as 
follows  (Figure  13)  with  numbers  right  justified. 


Data  Sheet 

Card  Colu 

1. 

1-5 

2. 

6-11 

3. 

12-16 

4. 

17-22 

5. 

28-30 

6. 

3  1-33 

7. 

34-36 

8. 

37-39 

9. 

(See 

Note  #1) 

40-41 

10. 

(See 

Note  #2) 

40-41 

11. 

42-44 

12. 

45-47 

13. 

48-50 
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14. 

(See  Note  #1) 

51-52 

15. 

(See  Note  #2) 

51-52 

16. 

(See  Note  #4) 

53-55 

17. 

(See  Note  #3) 

53-55 

18. 

56-58 

19. 

Does  not  apply 

20. 

(See  Note  #5) 

59-62 

21. 

(See  Note  #5) 

63-66 

22. 

71-72 

NOTES: 

1.  (a)  If  rhombic  is  circled,  enter  -1 

(b)  If  X/  2  dipole  is  circled,  enter-3 

(c)  If  vertical  is  circled,  enter  -2 

(d)  If  vertical  in  multiples  of  wave  length 
is  given,  enter  in  height  cols.  -14  for 
X/  4;  -12  for  X/  2;  -1  for  X,  etc. 

2.  If  this  number  is  given  instead  of  the 
information  which  appears  in  Note  1  (d) 
above,  enter  this  number. 

3.  (a)  If  check  in  Industrial,  enter  -1 

(b)  If  check  in  Residential,  enter  -2 

(c)  If  check  in  Rural,  enter  -3 

(d)  If  check  in  Remote  Unpopulous,  enter  -4 

4.  If  number  is  given,  enter  this  number  ini 
lieu  of  #17  from  data  sheet. 

5.  Bearing  in  degrees  east  of  north. 
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2.  Method  Card  Format 

The  following  is  a  list  of  methods  available.  The  numbers  on 
the  right  call  the  various  methods  when  placed  in  the  columns  shown. 


Card  Columns 

Calculations  to  be  Performed_ 2  4  6  8 

MUF-FOT  listing  from  circuits  on  cards  110  0 

MUF-FOT  listing  from  circuits  on  tape  10  0  0 


Sample 

Print- 

Out 

Figure  1 


FOT  curve  from  circuits  on  cards 
FOT  curve  from  circuits  on  tape 


1101  Figure  9 
10  0  1 


Reliability-mode -angle  from  circuits 
on  cards 

Reliability-mode -angle  from  circuits 
on  tape 

FOT-LUF  listing  every  hour  from 
circuits  on  cards 

FOT-LUF  listing  every  hour  from 
circuits  on  tape 


2  10  0  Figure  2 

2  0  0  0 


1  0  0  Figure  3 


FOT  listing  every  hour,  LUF  every  2 
hours  from  circuits  on  cards 
FOT  listing  every  hour,  LUF  every  2 
hours  from  circuits  on  tape 


3110  Figure  4 
3  0  10 


FOT-LUF  curves  from  circuits  on  cards 
FOT-LUF  curves  from  circuits  on  tape 


3  10  1  Figure  10 

3  0  0  1 


System  loss  from  circuits  on  cards 
System  loss  fro.m  circuits  on  tape 


4  1 

4  0 


0  0  Figure  5 

0  0 


Field  strength  (dbu)  from  circuits  on  4  1 

cards 

Field  strength  (dbu)  from  circuits  on  tape  4  0 


1  0  Figure  6 

1  0 


Reliability-FOT  every  2  hours  from 
circuits  on  cards 

Reliability-FOT  every  2  hours  from 
circuits  on  tape 


5  1 


0  0  Figure  8 


5  0  0  0 
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Available  S/  N  from  circuits  on  cards 
Available  S/N  from  circuits  on  tape 


6  10  0  Figure  7 

6  0  0  0 


3.  Assembling  the  Data  Decks 


Only  one  data  tape  is  needed  with  these  routines  (Tape  #2).  A 
utility  tape  #3  is  also  needed  which  will  contain  the  circuits  from  a 
given  run.  If  a  large  number  of  circuits  are  to  be  run  repeatably,  it 
is  wise  to  save  the  tape  and  run  from  it  in  lieu  of  circuit  cards. 

To  assemble  a  data  deck,  first  make  the  method  and  frequency 
complement  cards  for  the  desired  HF  predictions.  If  the  circuits  are 
to  be  read  from  cards,  the  circuit  cards  must  follow  the  frequency  card. 
Following  the  last  circuit  card,  insert  a  "nines  card"  (the  number  9 
punched  in  columns  1  through  19),  to  indicate  the  end  of  the  circuit 
cards.  The  next  cards  will  indicate  the  months  and  solar  activity 
levels  for  which  the  calculations  are  to  be  made. 

The  months  January  through  December  are  indicated  by  numbers 
1  through  12  placed  on  the  card  in  columns  23  and  24.  On  the  same 
card,  in  columns  25,  26  and  27,  punch  the  appropriate  sunspot  number. 
If  it  is  desired  to  terminate  the  calculations  of  these  circuits  under  the 
called  method,  a  card  with  -1  punched  in  columns  23  and  24  should 
follow  the  "month  and  sunspot  number"  card.  However,  if  another 
method  is  to  be  run  for  the  same  circuits,  follow  the  -1  card  by  method 
and  frequency  cards  that  allow  the  circuits  to  be  read  from  tape  (see 
method  card  format).  The  circuits  will  be  read  from  tape  #3  just 
prepared.  Following  the  new  frequency  and  method  cards,  place  the 
desired  month  and  sunspot  number  cards.  Notice  the  nine1  s  card  has 
been  omitted.  A  -1  card  inserted  after  the  month  and  sunspot  number 
card  will  terminate  the  method  and  a  blank  card  following  the  -1  card 
will  terminate  the  job. 
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The  data  deck  (to  be  run  from  cards)  is  composed  of  the  following 
cards  in  the  order  they  appear. 

1.  Method  card 

2.  Frequency  complement  card 

3.  Group  of  circuit  cards 

4.  Nines  card 

5.  Group  of  month  and  sunspot  cards 

6.  -1  card 

(See  Figures  14  and  15) 

Any  number  of  different  methods  with  different  circuits  may  be 
set  up  sequentially  in  the  above  order.  When  it  is  desired  to  terminate 
the  job,  place  a  blank  card  after  the  -1  card. 

A  data  deck  to  be  run  from  tape  is  composed  of  only: 

1.  Method  card 

2.  Frequency  complement  card 

3.  Group  of  month  and  sunspot  cards 

4.  -1  card 
(See  Figure  16) 


4.  Purpose  of  Cards  in  Data  Deck 


(a)  Method  Card: 

The  method  card  shown  in  Figure  13  allows  the  communicator  a 
choice  of  computations  shown  in  the  example  print-outs  (Figures  1 
through  10)  with  the  option  of  having  data  for  the  circuits  on  either 
cards  or  magnetic  tape. 

(b)  Frequency  Complement  Card: 

The  frequency  card  allows  the  communicator  a  choice  of  10  fre¬ 
quencies  for  which  calculations  will  be  made,  e.g.  ,  3.  1240,  4.  7856 
.  30.0000.  If  fewer  than  10  frequencies  are  wanted  merely  put 
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the  desired  ones  sequentially  on  card  from  the  beginning,  then  punch 
unused  fields  with  990000.  If  method  #3  is  to  be  run,  it  is  suggested 
that  a  frequency  card  identical  to  Figure  12  be  used. 

(c)  Circuit  Cards  (Figure  13): 

The  circuit  cards  are  designed  to  contain  all  necessary  circuit 
parameters  to  be  run  under  any  choice  of  the  computations  shown  in 
Figures  1  through  10. 

(d)  Nines  Card  (Figure  12): 

The  nines  card  is  a  circuit  card  which  is  not  computed.  The 
machine  senses  this  card  and  is  told  the  entire  list  of  circuits  has  been 
calculated  for  the  given  month. 

(e)  Month  and  Sunspot  Number  Card  (Figure  12): 

The  month  card  contains  the  month  of  the  year  and  its  associated 
sunspot  number  for  the  desired  circuit  calculations. 

(f)  Minus  One  Card  (Figure  12): 

The  minus  one  card  instructs  the  machine  that  all  computations 
are  complete  for  a  given  method. 

(g)  Blank  Card  (Figure  12): 

A  blank  card  following  a  "minus  one  card"  instructs  the  machine 
that  all  computations  for  a  given  run  are  complete. 


-  27  - 


1 

TRANSMITTER 
40.75N  -  73. 

93N 

JAN 

i 

32. 

SSN* 

RECEIVER 

33N  -  64.70VI 

25. 

135 

BR 

BEARINGS 
•8  321.3 

32. 

007 

N. MILES 
672.5 

6MT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

1 

5.7 

4.9 

7 

5.7 

4.9 

13 

10.3 

9.2 

19 

11.6 

10.8 

2 

5.2 

4.4 

e 

5.5 

4.7 

14 

11.4 

10.6 

20 

11.7 

9.9 

3 

5.2 

4.5 

9 

5.4 

4.6 

15 

11.6 

11.4 

21 

11.2 

9.5 

4 

5.4 

4.6 

10 

4.9 

4.2 

16 

11.8 

11.8 

22 

9.7 

8.2 

5 

5.4 

4.6 

11 

5.1 

4.3 

17 

11.8 

11.8 

23 

7.7 

6.6 

6 

5.7 

4.9 

12 

7.4 

6.9 

18 

11.5 

11.5 

24 

6.5 

5.5 

2 

TRANSMITTER 
51.50N  -  0. 

01N 

JAN 

1 

27. 

5SN* 

RECEIVER 

92N  -  15.67b 

25. 

212 

GC 

BEARINGS 
.5  22.2 

2. 

016 

N. MILES 
1581.3 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

1 

10.3 

8.8 

7 

8.1 

6.9 

13 

22.7 

19.3 

19 

13.5 

11.5 

2 

10.8 

9.2 

8 

12.9 

11.0 

14 

21.7 

18.5 

20 

11.8 

10.1 

3 

10.6 

9.0 

9 

19.6 

16.6 

15 

21.5 

18.2 

21 

10.8 

9.2 

4 

9.9 

8.4 

10 

22.8 

19.3 

16 

21.0 

17.8 

22 

9.8 

8.3 

5 

9.3 

7.9 

11 

23.6 

20.0 

17 

19.5 

16.6 

23 

9.7 

8.3 

6 

8.1 

6.9 

12 

23.7 

20.2 

18 

16.6 

14.2 

24 

10.1 

8.6 

3 

TRANSMITTER 
12.  ION  -  8. 

50E 

JAN 

1 

51. 

SSN- 

RECEIVER 

5 ON  -  0.01b 

25. 

351 

LO 

BEARINGS 
.8  167.0 

16. 

024 

N. NILES 
2400.4 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

GMT 

MUF 

FOT 

1 

13.4 

11.4 

7 

15.9 

13.6 

13 

25.0 

21.9 

19 

14.8 

12.6 

2 

13.0 

11.1 

8 

23.5 

19.9 

14 

25.3 

21.5 

20 

13.9 

11.8 

3 

11.8 

10.1 

9 

26.6 

22.6 

15 

24.4 

20.7 

21 

13.1 

11.1 

4 

10.7 

9.1 

1C 

27.3 

23.2 

16 

23.0 

19.5 

22 

13.0 

11.0 

5 

9.2 

7.8 

11 

27.4 

23.3 

17 

20.2 

17.1 

23 

13.0 

11.0 

6 

9.8 

8.4 

12 

26.6 

22.6 

18 

16.7 

14.2 

24 

13.0 

11.1 

Figure  1.  Computer  Print -Out  of  Circuit  MUF  and  FOT 
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Figure  8*  Alternative  Computer  Print-Out  of  Circuit  MUF,  FOT  and 
Circuit  Reliability 
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DATA  FOR  HIGH  FREQUENCY  COMMUNICATION  PREDICTIONS 


_  To  _  _ _ 

Transmitter  Location  Receiver  Location  Customer  Circuit  Number 


1.  Transmitter  Latitude  (Degrees  in  Decimal  Form)  North  =  +  ;  South  =  - 

2.  Transmitter  Longitude  (Degrees  in  Decimal  Form)  West  =  +  ;  East  =  - 

3.  Receiver  Latitude  (Degrees  in  Decimal  Form) 

4.  Receiver  Longitude  (Degrees  in  Decimal  Form) 

5.  Power  Delivered  to  Transmitting  Antenna 

Degrees 

Degrees 

Degrees 

Degrees 

KW 

6.  Transmitting  Antenna  -  Height 

7.  -  Leg  Length 

8.  -Tilt  Angle 

9.  -  Type;  Rhombic,  \/  2  Dipole,  or  Vertical 

10.  Transmitting  Antenna  Gain  if  Type  Other  Than  Item  #9 

Meter  s 

Meter  s 
Degrees 

DB 

11.  Receiving  Antenna  -  Height 

12.  -  Leg  Length 

13.  -  Tilt  Angle 

14.  -  Type;  Rhombic,  X./  2  Dipole,  or  Vertical 

15.  Receiving  Antenna  Gain  if  Type  Other  Than  Item  #14 

Meters 

Meter  s 
Degrees 

DB 

16.  Man-Made  Noise  at  Receiver  Input 

Frequency  Bandwidth  DB  <  1  Watt 

17.  Type  of  Receiving  Area  if  #16  is  Unknown 

Industrial  Residential  Rural  Remote  Unpopulous 

18.  Required  Hourly  Median  Signal-to-Noise  in  1  c/  s  Bandwidth 

DB 

19.  Detailed  Description  of  Service  if  #18  is  Unknown 

Type  of  Service,  e.g.  ,  Radiotelephone 

Modulation,  e.  g.  ,  SSB 

Bandwidth 

Number  of  Channels 

Words  Per  Minute  Per  Channel 

Maximum  Teletype  Error  Rate 

Minimum  Radiotelephone  Intelligibility 
_  Diversity  Employed 

Any  Other  Description  of  Service 

20.  Enter  Bearing  of  Transmitting  Rhombic  Antenna  if  Off -Great- 

Circle  Path 

21.  Enter  Bearing  of  Receiving  Rhombic  Antenna  if  Off-Great- 

Circle  Path 

22.  Two  Letter  Code  for  Receiver  Location 
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EXAMPLE  DATA  DECK  FOB  BUN  UNDER  TWO  METHODS  FROM  CARD 


II  -1 

III  6  34 

|99999|9|9999999999999 
|  1538^12052  2295-12020 
|  3S9?|  |  7468  3232  6473 
|  300pp  50000  70000 

I  3  1  111  0 


j 

10  32250  70-1  30180  67-1  -3  48  TM 

1  15  -2  66  96  67-1  -2  61  315  KU 

90000  110000  130000  150000  200000  150000  300000 


-1 

12  53 


11  60 

p9999pp999999999999 
|  1536J-I1 2050  3401  -7165 
|  3140]  f-7200  1520-12052 
|  3401]  |-7165  1530  -3900 


20  15  -3  25  96  64-1  -3  51 

40  25180  67-1  36103  66-1  -4  49  102  25 

10  21201  70-1-12  -2  -3  47 


'  I 


i  I 


5  1 
c 
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CKj 
AAj 


300pp  50000  70000  90000  100000  120000  150000  170000  180000  250000 


FORTRAN  STATEMENT 


n  i 

ill  1,  l 

i 

212  1 2  2 
I 

3i3  3  3  3 

i 

414  3  14 

I 

5'j  5  3  5 
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is  m  il  u  i4  is  ii  i;  ie  is  r  :  n  :i  it  n  3.  r  u  k  r  n  -v .a  -  * ;  it  k« <1  s;  4;  u  * «e  r  x  >  n  r  n  /  -  ■.  </ s  .  is n  K  c;  u  a  it  n  v]  r,  n  r  ;<  -  k 

1  I  I  I  I  I  I  1  1  1  I  1  i  I  1  II!  !  I  I  I  1  1  !  I  I  1  I  I  1  1  I  I  11  tl  I  !  1  II  11  1  1  i  1  I  1  1  1  4  (1  1  1  1  1  1  1  I  1  1  !il  1  1!  I  II  1 


OHRI  5T,:T|  T1TJTTSTTI  jWimTSTi  i  Ttlft  a  Cl  0  F  T  O  SI  D  Off  ffTFl  5  0  Cl  cT5  oil  c  i 
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3I3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  ?  3  3  3  jjs  3  3  3  3  3  3  3 
j 

4;4  444444  4  444444444444444444444444444444444444444444444444444444444  4j4  444444  4 
5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  S  5  5  5  5  5  5  5  5  5  5  5  5  5  E  5  5  5  5  5  5  5  E  5  5  5  5  5  5  5  5  5  E  5  5  5  5  E  2  5  sjj  E  5  E  E  5  E  5, 
$  6  6  6  6  6  C  $  6  6  6  C  6  6  6  S  5  E  B  6  6  6  6  C  6  6  5  c  i  5  6  5  E  6  6  5  E  6  C  E  6  6  E  S  C  E  6  E  E  E  E  C  B  6  E  E  6  C  8  E  E  E  E  6  E  sic  £  l  E  E  E  6  cj 
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8ifl  3  3  5  8:8  8  8  S  !  8  8  8  8  8  8  8  8  8  6  8  8  6  3  8  8  8  8  8  8  8  8  8  t  8  l  3  S  $  3  5  3  £  l  8  8  8  £  6  7  5  C  £  6  E  8  i  l  f  i  8  E  f  £  2  E  f  E  t  r  fc.G  t  f  £  fc  £  I  f! 


F igure  1 5 


|B-  3705 


-  42 


EXAMPLE  DATA  DECK  FOR  ONE  METHOD 
USING  CIRCUITS  FROM  MAGNETIC  TAPE 
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i 
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i 

2;2  2  2  2 
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I 
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FORTRAN  STATEMENT 


i 

3||  3  3 

I 

4i44  4  4 
I 

515  5  5  5 

r 

8  6  6  6  6 

I 

717  7  7  7 

I 


0 1  '0  ODD U  HOTO  D  DWdTiTOffDTOO  0  OTTTOTffTffTiTiroTirrD  0  0  A  0  (T£T 0  U  0  OTD  0  0  0  0  0  0  O'O  06  0  0  0  0  0  O'if 


l 

ii  i  io  n  u  n  it  u  is  u  ii  u  n  ii  a  n  ii  a  a  n  » is » ii  ini  it  s  »  «  a  » to  ti «  « tt  ts  ti  ti  ti  ti  eo  si  M  ei  it  a  s«  si  m  m  to  si  ei  m  ot «  « ei  h  ei  io  n  k 

n  i  n  u  n  n  i  ii  n  n  m  1 1  n  n  m  i  n  1 1  ii  1 1  ii  n  m  1 1  n  i  ii  n  i  m  i  m  i  ii  n 
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EXAMPLE  HF  PREDICTION  DECK  TO  EXECUTE  UNDER  IBM-7090  FORTRAN  MONITOR 


r.  r .... 


'  I  ft D <85330? MUFLUF. LUCAS. P> 
"i*  I  fr'RTr  1.0  'T5''46? 


^  $*ni  r-  s 

V;  •;»  i  .  % 

■ 


t4.f  ^l-iSWSS 


FORTRAN  STATEMENT 


T3T ITTl  0  0  ?  3  0  0  0  0  0  5  0  0  B  0  0  3  0  0  0  0  0  a  0  0  0  0  0  0  0  0  0  0  o  o oocoooooooooaoooooooooooooooo 

1 12  1  4  S  U!  7  J  9  IB  II  17  II 14  IS  IB  1 '  IB  19  TO  71  72  Z)  2«  75  75  77  78  73  7C  71  17  O  5«  ISIS  77  a  76 « 61  BIBS  t4 45  <6  67  *8  45 S3  SI  S7  57  MSS  MS?  H  IS  MSI  67  6)  67  65  66  67  68  63  10  71  77  71  76  75  76  77  a  78  81 
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i  1  I  i 
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1  1 1 
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Figure  17 
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VII.  PROGRAM  OUTPUT 

The  output  is  shown  in  Figures  1  through  10.  The  output  is 
designed  for  9  by  11  standard  computer  paper. 

A  description  of  variables  appearing  in  the  print-outs  reading 
left  to  right: 

1.  Sequential  circuit  number  appearing  in  the  extreme  upper 
left  is  generated  within  the  machine  and  is  designed  to  be  useful  for 
ready  reference  to  a  master  list.  Sample:  118 

2.  Month  of  year  for  which  predictions  are  calculated. 

Sample:  MARCH 

3.  Sunspot  number  used  in  predictions.  Sampe:  SSN  -  10. 

4.  Customer  identification  number  which  is  generated  from  a 
two-letter  code  for  receiver  location,  bearing  (receiver  to  transmitter), 
and  great  circle  distance  of  path.  Sample:  KD  12.  012 

5.  Transmitter  geographic  coordinates  in  hundredths  of 
degrees.  Sample:  13.45N  -  144.  80E 

6.  Receiver  geographic  coordinates  in  hundredths  of  degrees. 
Sample:  26.  72N  -  127.  80E 

7.  Great  circle  bearing  of  transmitter  to  receiver  and  bearing 
of  receiver  to  transmitter.  Sample:  3  12.4  •  126.5 

8.  Type  of  antenna  used  at  transmitting  terminal.  Sample: 
Rhombic  12H  39L  58  Deg. 

H  -  height  in  meters 
L  -  leg  length  in  meters 
Deg  =  tilt  angle  in  degrees 

9.  Level  of  man-made  noise  at  receiver  site.  Sample: 

NOISE  =  3 

10.  Type  of  antenna  used  at  receiving  terminal.  Sample: 
VERTICAL  -  14H  OL  ODEG. 


45  - 


The  -14H  indicates  the  grounded  vertical  is  1/4  (X).  The  same 
analogy  would  apply  for  -12H  (X/  2),  -1H  (X),  etc.  A  print-out  without 
the  minus  indicates  the  vertical  height  in  meters  (e.g.  ,  13H  indicates 
the  grounded  vertical  is  13  meters  high). 

11.  Great  circle  distance  of  path  in  nautical  miles.  Sample: 

N.  Miles  1243.2 

12.  The  effective  power  delivered  to  the  transmitting  antenna. 
Sample:  PWR  =  20.  00  kw 

13.  The  required  signal-to -noise  in  a  one -cycle  band  (db  >  1  watt) 
at  receiving  antenna  terminals.  Sample:  REQ.  S/N  =  45  db 

14.  Greenwich  mean  time.  Sample:  GMT 

15.  Maximum  Usable  Frequency.  Sample:  MUF 

16.  Optimum  Traffic  Frequency.  Sample:  FOT 

17.  Truncated  numbers  representing  the  frequency  complement 
card  (Figure  14),  Sample  4,  6,  8  .  .  .  .  23 

18.  Probable  mode  associated  with  least  loss.  Sample:  3F,  4E 
and  4X.  4X  indicates  propagation  by  E  and  F^  region  using  four  hops. 

19.  The  theoretical  angle  of  arrival  associated  with  the  most 
probable  mode  of  propagation.  Sample:  20 

20.  Reliability  defined  as  the  per  cent  of  days  with  the  month 
the  available  signal-to-noise  will  be  adequate.  Sample:  Reliability 

21.  Signal-to-noise  defined  as  the  available  monthly  median. 
Sample:  S/N  .  .  .  .  db 

22.  System  loss  defined  as  the  ratio  of  the  power  delivered  to 
the  transmitting  antenna  relative  to  that  available  at  the  receiving 
antenna  terminals.  Transmission  loss  may  be  obtained  by  using 
isotropic  antennas.  Sample:  Loss  .  .  .  db 

23.  The  Lowest  Useful  High  Frequency  defined  as  the  lowest 
frequency  that  will  give  90%  reliability  at  a  given  hour.  Sample:  LUF 
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24.  A  zero  appearing  in  the  LUF  column  denotes  a  LUF  greater 
than  the  FOT. 

2  5.  A  -3  in  the  LUF  column  denotes  a  LUF  below  3  Me/  s. 

26.  A  calculation  in  any  method  showing  all  zeros  denotes  a 
frequency  for  which  the  calculation  was  not  completed  due  to  mode 
limitations  or  the  primary  calculation  for  the  desired  mode  exceeded 
the  absolute  value  999. 
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VIII.  COMPLETE  FORTRAN  LISTING  OF  COMPUTER  PROGRAM 
FOR  PREDICTING  HF  SYSTEM  PERFORMANCE 

C  MUFLUF •  •  • *H— F  PREDICTIONS  BY  D*L. LUCAS  AND  G.W.HAYDON  »• •• • .CRPL. .« .. 

C  FORMATS  AND  STORAGE  ASSIGNMENTS 

6  FORMAT ( 10F7.4)  0001 

1  FORMAT (2(F5»2,F6*2) ,  I  2 , 5F3 .0 » 1 2 , 3F3 * 0 » 12 ♦ 1 3 »F3 *0 » 2F4* 0  * 4X , A2 )  0002 

3  FORMAT (10X*I5»10X»A6»  10X ♦ 4HSSN=F5 *  0 , 10X , A2 » F7 , 3/9X #  11HTR  0003 

2ANSMITTER»13X  » 8HRECE I VER  *12X*8HBEARINGS*6X*7HN#MILES/7X  ,F6.2,A1*  0004 

32H  -  »F7#2»A1,5X,F6*2,A1»2H  -  . F7 . 2 # A 1 #3X # 2F6. 1 #4X , F8 • 1 )  0005 

4  FORMAT  (F10.0)  0006 

5  FORMAT (412)  0007 

676  FORMAT ( 12A6 )  0008 

7  FORMAT ( 1X/6X  »  0009 

43HGMT»2X»3HMUF*3X»3HF0T»3(5X»3HGMT*2X ♦ 3HMUF  .3X.3HFOT)/  1H04(4X*I3*  0010 

52F6.1) /1H04(4X*I3*2F6.1 ) /1H04(4X* I3.2F6.1) /1H04(4X»I3»2F6*1)/1H04(  0011 

64X»I3*2F6«1)/1H04(4X*I3*2F6*1)/1H  Al)  0012 

DIMENSION  AMONI 12) »GLAT ( 5 ) ♦ ABI ( 24 ) » AB I Y I  5 » 24 ) »CLCK (5*24) »GY ( 5 ) •  0013 

lF2S(5f 24) »F2H(5*24) »GML(5»24) *GMH( 5*24)«MF(3) »ME(2)  0014 

2#ELD(10),FLD(10) » F24 ( 5 ) *GMA ( 5 ) » AD J < 4 ) » ABC (12)  0015 

DIMENSION  RD(5)»SUN(12)*A(10»7»14)»RASSN(12)*EMF(5*24)  0016 

DIMENSION  0( 20»60*4) , IL(4) »JL(4) »KL(4) »LK(4) »JAL(4)  »CKC(24)  0017 

1»G(60) ,AB( 60) #5(20,24) »C< 20 »24) .GAMMA ( 5 , 24 ) ,3A ( 60 )  0018 

DIMENSION  CLAT(5) ,CLONG(5) ,EMUFY(5,24) #FMUFY(5#24),  UFY ( 5 ) #FOTY ( 5 )  0019 

5  #  I G ( 24 )  0020 

DIMENSION  UF(24) ,F0T(24) »P(29.16,6) #ABP( 2#6)  0021 

DIMENSION  AHA(2) » AR A ( 2 ) # AVA ( 2 )  0022 

DIMENSION  XI XX ( 1 2  )  0023 

DIMENSION  FREL ( 12 )  0024 

COMMON  ABP.P  0025 

COMMON  A,0,G»AB#S,C»GAMMA,BA,EMUFY,FMUFY#IG»UF#F0T,SUN#AVA#ARA,GLA  0026 

1T.ABI »AB!Y,CLCK»GY»F2S»F2H,GML,GMH,MF,ME,ELD#FLD.F24,GMA,ADJ,ABC#  0027 

2EMF#CKC,CLAT,CL0NG,AK,BK,CK.DK,PI2»EK,EEK,GLT,GLG,AAA,ALA»ASA,  0028 

3AFCiANC#AWC.AEC#ASC#AZZ#MAN#RSN#<W*PWR#IANT»IANR,Y2»X2#MOUSE»SSN#  0029 

4N0CIR*K4,  AX.XLONG.AL ,YLAT»XLAT,YL0NG»YL»BTRY,GCD<M,GCDNM#K1,K2#K  0030 

53»<5#GCD,IRSN,M0N,ID»M0USE#AY,METH0D#AHA,BM0NS»MAP  0031 

COMMON  FREL  0032 

COMMON  XNH,XNL#XND,RNh,RNL,RND»HA»JIG»HAR#XTR,XETA,WETA»IHR#MIT  0033 

REWIND  2  0034 

REWIND  3  0035 

909  FORMAT ( 1H1 )  0036 

DO  101  !«1»24  0037 

101  IGID-I  0038 
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C  CONSTANTS  TO  BE  USED  FREQUENTLY 


AK  =  1 .745329E-2 

0039 

BK=  5 • 729577E 1 

0040 

CK=1 *  1 1 136E2 

0041 

DK= • 0062137 

0042 

P I 2* 1 *57079 

0043 

EK=. 539956 

0044 

EEK=.62i4 

0045 

GLT=78.*AK 

0046 

GLG=70.*AK 

0047 

B 

AMON< 1 >=412145606060 

0048 

B 

AMON ( 2 ) =262522606060 

0049 

B 

AMON(3 >=442151606060 

0050 

B 

AMON (4 >=214751606060 

0051 

B 

AMON ( 5  >  =442170606060 

0052 

B 

AMON ( 6 ) =416445606060 

0053 

B 

AMON ( 7  >  *416443606060 

0054 

B 

AMON ( 8  > =216427606060 

0055 

B 

AMON ( 9 ) =622547606060 

0056 

B 

AMON  ( 10=462363606060 

0057 

B 

AMON (11) =454665606060 

0058 

B 

AMON ( 121=242523606060 

0059 

B 

AAA= 60 606060 60 60 

0060 

B 

AFC *26606060 6060 

0061 

B 

AHA ( 1 ) =606060603040 

0062 

B 

AHA( 2 >=243147464325 

0063 

B 

AVA< 1) =606060606525 

0064 

B 

AVA( 2 >=516331232143 

0065 

B 

ARA ( 1 ) *606060606051 

0066 

B 

ARA ( 2 ) *304644223123 

0067 

B 

AZZ= 00 6060 6 06060 

0068 

B 

ANC =456060606060 

0069 

B 

AWC= 66 6060 606060 

0070 

6 

AEC= 256060606060 

0071 

B 

A SC =6 26 0606 060 60 

0072 

B 

FNINE=111111111111 

0073 

SUN ( 1 ) =-21 • 

0074 

SUN(2)=-13. 

0075 

SUN ( 3 ) *~03 • 

0076 

SUN ( 4 ) =09» 

0077 
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SUN ( 5 ) *  18  #  0078 

SUN ( 6 ) =23*  0079 

SUN ( 7 ) =21 •  0080 

SUN(8)»14.  0081 

SUN { 9 ) =04*  0082 

SUN ( 10 ) =-08 •  0083 

SUN ( 11 ) =— 1 4 •  0084 

SUN ( 12 ) =-23*  0085 

ADJ ( 1 ) =9*  0086 

ADJ ( 2 ) = 17 •  0087 

ADJ (3) =20.  0088 

ADJ ( 4 ) =28 •  0089 

C  READ  CIRCUIT  DATA  AND  SEARCH  DATA  TAPE 

WRITE  OUTPUT  TAPE  6.909  0090 

464  READ  INPUT  TAPE  5.5,  0091 

1  METHOD. METH. IHR. MAP  0092 

READ  INPUT  TAPE  5*6. ( FREL ( I ) . 1=1 . 10 )  0093 

M I T= 1  0094 

IF  (METHOD)  463.12.463  0095 

12  REWIND  2  0096 

WRITE  OUTPUT  TAPE  6.909  0097 

REWIND  3  0098 

CALL  EXIT  0099 

463  L 1=2  0100 

GOTO  (30.30.30.30.31.30) .METHOD  0101 

30  L 1=1  0102 

31  NOC I R=  0  0103 

MONTH=0  0104 

IF(METH)  7770.7771.7770  0105 

7770  READ  INPUT  TAPE  5.676.X1XX  0106 

WRITE  OUTPUT  TAPE  3.676»X1XX  0107 

IF(FNINE-X1XX( 1 ) )  7770,7771.7770  0108 

7771  REWIND  3  0109 

85  READ  INPUT  TAPE  5.1.  0110 

1  X1X,Y1Y.X2X,Y2Y,MONS.SUNS.TWR,XNhO,XNLO,XNDO. ITANT.RNH0.RN  oill 

1LO.RNDO. IRANT.MANN.RRSN.XETA.WETA.HAO  0112 

10  X1=X1X  0113 

Y1=Y1Y  0114 

X2=X2X  0115 

Y2=Y2 Y  0116 
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I  ANT  =  I TANT  0117 

I  ANR  = I  RANT  0118 

MAN=MANN  0119 

RSN=RRSN  0120 

I RSN=RSN  0121 

PWR*TWR  0122 

XNH=XNHO  0123 

XNL=XNLO  0124 

XND=XNDO  0125 

RNH=RNHO  0126 

RNL  =RNLO  0127 

RND=RNDO  0128 

HA=HAO  0129 

IF(MONS)  464  .11,82  0130 

82  LOCK=  MONS-MONTH-1  0131 

SSN=SUNS  0132 

MONTH=MONS  0133 

MOUSE=MONS  0134 

KZ=1  0135 

READ  INPUT  TAPE3*l*XlX.YlY.X2X.Y2V»MONStSUNStTWR»XNHO.XNLO»XNDO»  0136 

1 ITANT.RNHO.RNLO.RNDO. I  RANT .MANN.RRSN.XETA.WETA.HAO  0137 

GO  TO  10  0138 

11  GO  TO  (  15.77  )  ,  KZ  0139 

15  NOC I R=Q  0140 

WRITE  OUTPUT  TAPE  6.909  0141 

62  J IG=0  0142 

M I T= 1  0143 

KZ=2  0144 

I F ( LOCK )  77,731.33  0145 

33  D0730  I U= 1 . LOCK  0146 

730  READ  TAPE  2  0147 

731  READ  TAPE  2. IL.JL.KL.LK.JAL.Q.A.P.ABP  0148 

77  READ  INPUT  T APE3 . 1 »X1 X . Y 1 Y ,X2X , Y2 Y . MONS . SUNS , T WR . XNHO .XNLO ,XNDO .  0149 

1 ITANT.RNHO.RNLO.RNDO. IRANT.MANN.RRSN.XETA.WETA.HAO  0150 

IFI91.-X1X)  78,78.16  0151 

78  REWIND  3  0152 

READ  INPUT  TAPE  5.1.  0153 

1  X  1  X.Y1Y.X2X.Y2Y. MONS .SUNS. TWR.XNHO.XNLO.XNDO, I T ANT . RNHO , RN  0154 

1LO.RNDO, IRANT.MANN.RRSN.XETA.WETA.HAO  0155 

I F ( MONS )  83.16,83  0156 
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83  I F ( MON S -MONTH )  84,16,15  0157 

84  REWIND  2  0158 

M0NTH=0  0159 

16  SSP=SUN ( MOUSE ) *AK  0160 

BMONS=AMON (MOUSE )  0161 

NOC I R=NOC I  R  +  l  0162 

JIP=1  0163 

JEP= 1  0164 

C  GREAT  CIRCLE  DISTANCE  AND  DISTANCE  TO  CONTROL  POINTS 

XX 1 =X 1 *AK  0165 

XX2=X2*AK  0166 

YY1=Y1*AK  0167 

YY2=Y2*AK  0168 

GCD=SINF (XXI ) *S I NF ( XX2 )+COSF( XXI )*COSF ( XX 2 ) *COSF( YY2-YY1 )  0169 

GCD=ACOSF ( GCD )  0170 

RD ( 1 >  =GCD/ 2 •  0171 

RD(2)=18.*AK  0172 

RD(3)=GCD-18.*AK  0173 

RD(4)=9.*AK  0174 

RD(5)=GCD-9.*AK.  0175 

GCDKM=GCD*C<*BK  0176 

GCDNM=GCD<M*EK  0177 

C  BEAR  I NGSS —  FORWARD  AND  BACKWARD 

DF=YY2~YY1  0178 

IF(ABSF(DF)-3. 141592654)  26,22,22  0179 

22  IF(DF)  24,25,25  0180 

24  DF=6.283185308+DF  0181 

GO  TO  26  0182 

25  DF=-6» 283185 308+DF  0183 

26  U=.5*( PI2-XX1+PI 2-XX2+GCD )  0184 

CIND=SINF(U)*SINF(U-PI2+XX1)  0185 

IF(CIND)  27,28,27  0186 

27  BTRY=114.5816*ATANF(SQRTF(ABSF(SINF(U-PI2+XX2)*SINF(U-GCD) ) / C I ND ) )  0187 

IF(DF)  7029,21,39  0188 

21  IFIXX1-XX2)  28,39,39  0189 

28  BTRY=180.  0190 

GO  TO  39  0191 

7029  BTRY=360.-BTRY  0192 

39  RF= YY 1-YY2  0193 

IF(ABSF(RF)-3. 141592654)  126,122,122  0194 
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122  IF(RF)  124.125,125  0195 

124  RF=6.283185308+RF  0196 

GO  TO  126  0197 

125  RF=-6.283185308+RF  0198 

126  U=.5*( PI2-XX2+PI 2-XX1+GCD)  0199 

CIND=SINF(U)*SINF(U-PI2+XX2>  0200 

IF(CIND)  57,778.57  0201 

57  XTR=114.5816*ATANF(SQRTF< ARSF < S I NF < U-P I 2+XX 1 ) *S I NF ( U-GCD ) )/CIND>  >  0202 

IF(RF)  8029,51,79  0203 

51  IFIXX2-XX1 )778,79,79  0204 

778  XTR=180.  0205 

GO  TO  79  0206 

8029  XTR  =360. -XTR  0207 

79  I F ( 4000 . -GCDKM )  40,41,41  0208 

40  <1  =  1  0209 

<2  =  2  0210 

K  3  =  3  0211 

K4=4  0212 

<5*5  0213 

GO  TO  43  0214 

41  <1=1  0215 

<2*1  0216 

<3*1  0217 

<4*1  0218 

<5*1  0219 

C  CONTROL  POINT  LOCAL  TIME  .  GEOGRAPHIC  LATITUDE  AND  LONGITUDE 

43  DO  61  IT=L1,24,L1  0220 

GMT  =  IT  0221 

COLO  =  GMT-Y  2/ 15 .  0222 

I F ( 24.— COLO )  94,95,95  0223 

94  COLO* COLO- 2 4 •  0224 

GO  TO  96  0225 

95  I F (COLO )  98,98,96  0226 

98  COLO=COLO+24 •  0227 

96  C<C< I T ) =COLO  0228 

SSL* ( 15.*GMT-180. )  0229 

DO  61  I  I =<1 *<5  0230 

PP=RD ( I  I )  0231 

CENLAT=COSF( PP)*COSF(PI2~XX2 )+SINF(PP )  0232 

l*SINF(PI2-XX2)*COSF(BTRY*A<)  0233 
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CENLAT=(PI2-AC0SF(CENLAT) )*BK  0234 

CENLG=  ( (COSF ( PP )-COSF ( PI 2-XX2 ) *C0SF ( P I2-CENLAT*AK ) )/(SINF(PI2  0235 

1-XX2  >*SINF(PI2-CENLAT*AK)  ) )  02  36 

IFIABSF (CENLG)-l, )  3001,3001,3000  0237 

3000  CENLG= 1 •  0238 

3001  CENLG  =  ACOSF { CENLG  )  0239 

CENLG= ( YY2-SIGNF (CENLG,DF ) )*B<  0240 

I F ( 180,-ABSF (CENLG) )  71,69*69  0241 

71  I F (CENLG )  73,72,72  0242 

73  CENLG=360.+CENLG  0243 

GO  TO  69  0244 

72  CENLG=CENLG-360.  0245 

C  GEOMAGNETIC  LOCATION  OF  CONTROL  POINTS 

69  GAT=ACOSF(SINF(GLT )  *S  I  NF  (  CENLAT*A<  ) +COSF  (  GLT  )  *COSF  <  CENLAT*AK.  )  *  0246 

1C0SF(CENLG*AK.-GLG)  )  0247 

GLAT ( II  ) =P I  2 -GAT  0248 

CLAT (II)  =CENLAT  0249 

70  CLONG ( I  I ) =CENLG  0250 

C  SUNS  ZENITH  ANGLE 

Z=<SSL-CENLG)*AK  0251 

56  CYCEN  =  SINF(CENLAT*AK  )*SI  NF  (  SSP  )+COSF  (  CENLAT*AK. )  *COSF(  SSP  )  *COSF  (Z  )  0252 

CYCEN=ACOSF(CYCEN)  0253 

CYCEN=ABSF(CYCEN)  0254 

C  ABSORPTION  INDEX  I  AT  CONTROL  POINTS 

32  I F ( 102, *AK-C YCEN )  58,58,59  0255 

58  AB I Y (  I  I  ♦  IT  )  =0,  0256 

GO  TO  886  0257 

59  ABIY( I  I , IT )  =  ( 1 .+  .0037*SSN ) *COSF ( .881*CYCEN ) **1 . 3  0258 

886  GO  TO  (165, 166), K2  0259 

165  ABI ( IT ) =ABIY( I  I , IT )  0260 

GO  TO  61  0261 

166  ABI (  IT  )  =  ( ABI Y( 1 , IT J+ABI Y( 2* IT )+ABI Y( 3, IT )  ) /3.  0262 

61  CONTINUE  0263 

DO  2112  II =K 1 , K5  0264 

DO  2112  IT=L1,24,L1  0265 

ABIC  =  AB I Y ( I  I ♦ I T )  0266 

C  E-LAYER  DISTANCE  FACTOR 

ARC  =  DIC*GCDKM  0267 

I F ( 16, -ARC )  88,87,87  0268 

87  E LFC  =  (  (  < ( (  ( -4, 368460907 E“9*ARC+ 1,3 344942 6 1 E-7 )*ARC-5. 9766 18436E-6 )  0269 
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1 *ARC+2. 6248083 1 5 E~4 ) *ARC~5 .038476266E-3 ) *ARC+3 . 761 38505 3 E-2 ) *ARC-1  0270 

2.133200756E-2)*ARC+.2085  0271 

GO  TO  100  0272 

88  ELFC= 1 • 02  0273 

C  E-2000  CONTROL  POINT  MUF 

100  EMC= ( ( <-l0.66484998*ABIC+39.26151056)*ABIC-52.41191754)*ABIC+37.67  0274 

1726072 )*ABIC+3. 345996232  0275 

EMF ( 1 1 » I T ) =EMC  0276 

2112  EMUFY (II*IT)=EMC*ELFC  0277 

C  GENERATE  CONTROL  POINT  FOF2  AND  M-3000  FACTORS 

DO  1112  1 1  =  K2  »K3  0278 

504  CENLG=CLONG (II)  0279 

CENLAT =CLAT (II)  0280 

IF(CENLG)  412.414.414  0281 

414  CLG=360.-CENLG  0282 

GO  TO  413  0283 

412  CLG=ABSF(CENLG)  0284 

413  BOY=SINF (CL G * A K )  0285 

COG=SINF(2.*CLG*AK)  0286 

DOG=COSF ( 2.*CLG*AK )  0287 

GOB=COSF(CENLAT*AK )  0288 

BOG=SINF(CENLAT*AK )  0289 

HOG=COSF(CLG*AK)  0290 

DO  500  10=1.4  0291 

415  I = I L ( 10)  0292 

J=JL (10)  0293 

K=KL (10)  0294 

L  =  LK ( 10)  0295 

JA= JAL ( 10 )  0296 

DO  408  ICAal.K  0297 

KP=KA-1  0298 

408  G(KA)=BOG**KP  0299 

LA=0  0300 

LO=L-l  0301 

KK=K+1  0302 

DO  409  KA  =  KX  » LO. 2  0303 

G ( KA ) =BOG**LA*GOB*HOG  0304 

GIKA+l )=BOG**LA*GOB*BOY  0305 

409  LA=LA+ 1  0306 

LB=0  0307 
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409  LA=LA+ 1  0306 

LB=0  0307 

LL=L+1  0308 

IM=I-1  0309 

DO  410  KA  =  LL  » I M  *  2  0310 

G( KA ) =BOG**LB*GOB**2*DOG  0311 

G(KA+1 ) =BOG**LB*GOB**2*COG  0312 

410  LB=LB+1  0313 

DO  411  JB=1 * J  0314 

AB ( JB ) =0 ,  0315 

I S=2* JB-1  0316 

DO  411  KA=1»I  0317 

411  AB( JB)=AB( JB)+Q< IS , KA , 1 0 ) *G ( KA )  0318 

DO  407  JB  =  2  » J  0319 

B A ( JB ) =0  «  0320 

IS=2*JB-2  0321 

DO  407  KA=1 » I  0322 

407  BA( JB) =BA( JB)+Q(IS,KA, IO)*G(KA)  0323 

DO  500  IT=L1,24*L1  0324 

503  GMT= IT  0325 

CLOCK=GMT-CENLG/15.  0326 

I F ( 24 • -CLOCK )  64,65*65  0327 

64  CLOCK  =  CLOCK-24 *  0328 

GO  TO  66  0329 

65  I F ( CLOCK )  68,68,66  0330 

68  CL0CK=CL0CK+24,  0331 

66  CLCK  ( II »IT)=CLOCK  0332 

TIME=(15.*CLOCK-180. )*AK  0333 

DO  998  JB  =  2  » J  0334 

FB= JB- 1  0335 

S( JB»IT)=SINF(FB*TIME)  0336 

998  C( JB,IT)=COSF(FB*TIME)  0337 

GAMMA( IO,IT)=AB(l)  0338 

DO  500  JB  =  2  »  J  0339 

500  GAMMA ( 1 0 , I T ) =GAMMA ( 1 0 , I T ) +AB ( JB ) *C ( JR ♦ I T ) +BA ( JB ) *S ( JB , I T )  0340 

F2-REGION  GYRO-FREQUENCY 

132  CALL  POLY ( 1 ,7, 7*CEN LG/180, *CENLAT /90* ,Y)  0341 

GYRO=Y / 2 •  0342 

GY (II) =Y  0343 

F2-LAYER  CONTROL  POINT  F2-4000  MUF 
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F2S( 1 1 • IT) =F2LS  0346 

F2HS=GAMMA ( 2 . I T ) *GAMMA ( 4» I T ) *1 •  1  0347 

F2H< II »IT)=F2HS  0348 

C  F2-LAYER  DISTANCE  FACTOR 

IFI24.-ARC)  1 57  1 1 57  1 1 58  0349 

157  FLFC=1 •  0350 

GO  TO  159  0351 

158  FLFC= ( «  <<< (-6.712654756E-9*ARC+4.49151441E-7)*ARC-9.985831104E-6)*  0352 

lARC+6.865259817E-5)*ARC+9.202437  3  32E-5)*ARC-*-2.264634341E-3)*ARC+4,  0353 

2699243101E-3  )*ARC  0354 

159  GAMMA ( 1 » I T ) =GAMMA ( 1 » I T ) +GYRO  0355 

GAMMA  <  2  1 1  T  ) =GAMMA ( 2  » I T )+GYRO  0356 

FMLS  =  GAMMA (1»IT)+FLFC*( F2LS-GAMMA ( It  IT) )  0357 

GMU  1 1  t  IT)  *GAMMA(  1 1  IT  )  0358 

FMHS  =  GAMMA (2#IT)+FLFC*( F2HS-GAMMA ( 2  1 1 T ) )  0359 

GMHII I t IT) =GAMMA< 2t IT )  0360 

C  F2-LAYER  CONTROL  POINT  MUF 

FMUFY ( I  I 1 1 T )  =  ( FMLS*( 180.-SSN ) +FMHS* ( SSN-10 . ) ) /170.  0361 

1112  CONTINUE  0362 

C  CIRCUIT  MUF  AND  FOT 

DO  1012  IT=Llt24,Ll  0363 

GMT  = I T  0364 

DO  1000  I  I =K2  * K3  0365 

IFIK2-1)  602t603t602  0366 

603  I R= 1  0367 

GO  TO  604  0368 

602  I R= I  I +2  0369 

604  I F ( EMUF Y ( I R  *  I T ) -FMUFY (  I  I  t  I  T  )  )  1007 1 1007 1 1008  0370 

1007  UFY( II )«FMUFY( I I.IT)  0371 

GO  TO  1009  0372 

1008  UFY( I  I >=EMUFY(  IRt  IT)  0373 

1009  FOTFY=.85*FMUFYl I  I t IT )  0374 

IF( (EMUFYI IRtIT)-FOTFY) )  lOlOt 101 1 1 1011  0375 

1010  FOTY ( I  I ) =FOTFY  0376 

GO  TO  1000  0377 

1011  FOTYI I  I ) =EMUFY( IR 1 1 T )  0378 

1000  CONTINUE  0379 

1 1 =K2  0380 

IF(l-II)  606  t605  1 606  0381 

605  UF (IT) =UFY ( 1 )  0382 
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FOT  (  IT)  =FOTY ( 1 )  O303 

GO  TO  1012  0384 

606  IF  (  UFY(2)~  UFY(3))  607*607.608  0385 

607  UF ( I T ) =  UFY ( 2 )  0386 

GO  TO  609  0387 

608  UF ( I T ) =  UFY ( 3 )  0388 

609  IF  (FOTY(2)-FOTY(3)  )  610,610.611  0389 

610  FOT ( I T ) =FOT Y ( 2 )  0390 

GO  TO  1012  0391 

611  FOT ( I T ) =FOTY ( 3 )  0392 

1012  CONTINUE  0393 

IRY=BTRY/10.  0394 

TRY= I R Y  0395 

HARaTRY+.OOOO 1*GCDNM  0396 

XLAT  =ABSF ( XI )  0397 

XLONG* ABSF ( Y 1 )  0398 

YLAT  =ABSF ( X  2 )  0399 

YLONG=ABSF(Y2)  0400 

AX= ANC  0401 

I  F ( X 1 )  700,701,701  0402 

700  AX= ASC  0403 

701  AL=AWC  0404 

IF(Y1)  702,703,703  0405 

702  AL=AEC  0406 

703  AY=ANC  0407 

I F ( X2 )  704.705,705  0408 

704  AY  =  ASC  0409 

705  YL=AWC  0410 

I  F ( Y2 )  706.2221.2221  0411 

706  YL=AEC  0412 

2221  GO  TO  (2223.2224,2224,2224,2224,2224) .METHOD  0413 

2223  I F (MAP  )  3223,322  3.3224  0414 

3224  CALL  CURVY ( FOT ,UF )  0415 

GO  TO  10  0416 

3223  WRITE  OUTPUT  TAPE  6 , 3 . NOC I R , BMONS ,SSN , HA »HAR , XL AT , AX , XLONG , AL ,  0417 

1 YLAT , AY , YLONG ,YL,XTR,BTRY,GCDNM  0418 

WRITE  OUTPUT  TAPE  6,7,  0419 

1  ( I G ( I ) »UF( I ) »FOT ( I ) , IG< 1+6 ) »UF ( 1+6 ) »FOT ( 1+6 ) . IG( 1+12 ) ,  0420 

2  UF ( 1+12) »FOT( 1+12) »IG( 1+18) »UF( 1+16 ) ,FOT( 1+18 ) .1=1 .6) .AAA  0421 

I F ( 3-MI T )  969,969.970  0422 
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969 

WRITE  OUTPUT  TAPE  6,909 

0423 

MI  T=1 

0424 

GO  TO  10 

0425 

970 

MI  T  =  MI T  +  l 

0426 

GO  TO  10 

0427 

2224 

CALL  LUFFY 

0428 

GO  TO  10 

0429 

END 

0430 
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SUBROUTINE  LUFFY  0001 

FORMATS  AND  STORAGE  ASSIGNMENTS 

DIMENSION  FDEK{ 10) ,EOEK ( 10)  0002 

DIMENSION  AMON( 12 ) ,GLAT( 5) » AB I ( 24 ) ♦ AB I Y ( 5 » 24 ) »CLCK( 5*24) *GY( 5 ) »  0003 

1F2S(5»24)»F2H(5*24) ,GML( 5*24) »GMH( 5 .24) *MF(3) *ME( 2) ♦ ELL ( 10) »FLL(10  0004 

2)  »ELD( 10) »EDEL{ 10) »FLD( 10) »FDEL( 10) »F24( 5) »GMA ( 5 ) *FLF ( 10 ) *FSKY ( 10 )  0005 

3  #  ESKY ( 10) ,GE(10) ,GF( 10) ,FSE( 10) , FSF ( 10 ) *  NOS ( 11 ) *FXMT( 10) *EXMT( 10) ,  0006 

4FRCR (10) ,ERCR( 10 ) , ADJ (4)  »ABC (12  5  *MODF( 12 )  0007 

DIMENSION  Q( 20*60*4) *IL(4) »JL(4) *KL(4) ,LK<4) *JAL(4)  *CKC(24)  0008 

1 *G( 60) *AB( 60) *S<  20*24) *C( 20*24) *GAMMA( 5*24) *BA (60 )  0009 

DIMENSION  RD( 5) *SUN( 12 ) »A( 10»7, 14) ,RASSN( 12 ) »EMF( 5,24 )  0010 

DIMENSION  CLAT (5 ) »CLONG ( 5) »EMUFY( 5*24) » FMUFY ( 5*24) »  UFY ( 5 ) » FOTY ( 5 )  0011 

5*IG(24)  0012 

DIMENSION  UF(24) ,F0T(24) *FF(24) *P(29*16*6) »ABP(2.6)  0013 

DIMENSION  AHA(2) »ARA( 2 ) *AVA( 2 ) »RANT( 2 ) *TANT( 2)  0014 

DIMENS  I ONNXLOS (11) ♦NANGLE(ll)  0015 

DIMENSION  TOP ( 4 )  0016 

DIMENSION  IFR(IO)  0017 

DIMENSION  FR EL ( 1 2 )  0018 

COMMON  ABP  *P  0019 

COMMON  A*Q»G*AB»S»C»GAMMA»BA»EMUFY  *FMUFY  » IG»UF  »FOT  *SUN*AVA*ARA»GLA  0020 

1T»ABI »ABIY,CLCK»GY,F2S»F2H,GML»GMH,MFtME*ELD*FLD»F24,GMA,ADJ,ABC,  0021 

2EMF*CKC>CLAT»CL0NG,AK,BK,CK»DK,PI2»EK»EEK,GLT»GLG*AAA»ALA*ASA*  0022 

3AFC »  ANC  »  AWC »  AEC »  ASC » AZZ  *MAN  *RSN  » KW » PWR  » I  ANT ♦ I ANR  » Y2  * X2  » MOUSE  » SSN  *  00  23 

4N0CI R  » K4 »  AX  » XLONG  *  AL ♦ YL A T » XL AT » YLONG » YL , P TR Y  » GCD<M , GCDNM , K 1 , <2 ♦ K  0024 
53*K5*GCD» I RSN  »MON  » ID  » MOUSE  »AY*METHOD»AHA  »BMONS  *MAP  00  2  5 

COMMON  FREL  0026 

COMMON  XNH  »XNL  *XND*RNH  *RNL  »RND»  HA  »  J IG  *HAR  » XTR  *  XETA  * WET  A ♦ IHR  *  M IT  0027 

COMMON  KAP»KUP*MA,TANT» INH» INL* IND»RANT» JNH» JNL» JND  0028 

2  F0RMAT(2X*2A6#1X»I3»1HH»1X*I3»1HL»1X*I3»3HDEG»12X»6HN0ISE=I3»12X»4  0029 

1HANT  =  I  3  »  2HDB )  0030 

1  FORMAT ( 2X*2A6*1X» 1 3 , 1HH » 1 X » I  3 . 1HL ♦ 1 X , 1 3 » 3HDEG » 4X * 6HNO I SE  =  I  3 , 2 A6 , 1 X  00  31 
1 *I3*1HH.1X*I3*1HL*1X* I3»3HDEG)  0032 

80  FORMAT (6X»4HPWR=F6.2*2HKW»11X»21H0PERATING  FREQUENC I ES , 9X , 14HF I  ELD  0033 
1  STRENGTH/5X »14HGMT  MUF  FOT , IX ) 1 0 1 4 ,4H  FOT )  0034 

3  FORMAT ( 10X» I5*10X*A6*  10X , 4HSSN=P5 . 0 ♦ 10X » A2 ♦ F7 , 3 /9X »  11HTR  0035 

2ANSMITTER,13X»8HRECEIVER* 12X  * 8HBEAR I NGS » 6X , 7HN , MI LES /7X  *F6,2  *A1  *  0036 

32H  -  *F7.2,A1»5X,F6.2,A1,2H  -  ♦  F7 . 2 » A 1 , 3 X  ♦  2 F6 . 1 » 4X , F8 . 1 )  0037 

4  FORMAT (6X»4HANT=I3»2HDB»13X*  6HNO ISE=I3»7X»2A6»1X»I3»1HH»1X»I3*1HL»  00  38 

11X* I3»3HDEG)  0039 


-  60  - 


9  F0RMAT(4X» I3.2F6. 1/20X.11 < I3.A1 )  ,6H  MODE /20X ,  1 1 1  4 . 7H  ANGLE /4X  »  0040 

213.2 F6.1.1X.11I4.13H  RELIABILITY)  0041 

12  FORMAT! IX. A2.F7. 3*15, 2F10. 2)  0042 

13  FORMAT(6X*4HANT=I3»2HDB.20X.6HNOISE=I 3*20X,4HANT=I3.2HDB)  0043 

14  F0RMAT(6X»4HPWR=F6.2»2HKW.11X,21H0PERATING  FREQUENC I ES . 10X . 8HREQ .S  0044 

1 /N= I  3  » 2HDB/5X  » 14HGMT  MUF  FOT *  IX » 1 0  1 4  * 4H  FOT )  0045 

29  FORMAT (4X.I3.2F6.1/20X.11! I  3  *  A 1 ) »6H  MODE/20X  » 1 1 1  4  *  7H  ANGLE /4X  *  0046 

2I3*2F6*1 *  1 X  *  1 1 14.9H  S/N..DB)  0047 

34  FORMAT(6X»4HPWR  =  F6.2.2HKW.16X.13HRELIABILITIES*13X.8HREG).S/N=I3*2H  0048 

1DB/9X. 14HGMT  MUF  FOT ♦ 1 X » 10 1 4 »9H  FOT  MCS)  0049 

7  FORMAT  <6X.  0050 

43HGMT  »2X»3HLUF»3X»3HFOT  *3( 5X*  3HGMT  *2X*3HLUF»3X»3HFOT ) /  1H04I4X.  13*  0051 

52F6.1)/1H04(4X,I3»2F6.1)/1H04(4X, I3.2F6.1 )/lH04(4X. I3.2F6.1)/1H04(  0052 

64X.I3.2F6.1) /1H04(4X,I3.2F6.1 )/lH  Al)  0053 

8  FORMAT! 31X.21HOPERATING  FREQUENC I E S/ 5X , 1 4HGMT  MUF  FOT. IX. 1014*  0054 

14H  FOT)  0055 

19  F0RMAT(4X. I3.2F6.1/20X.11 ( I3.A1 ) .6H  MODE /20X . 1 1 1  4 » 7H  ANGLE /4X  .  0056 

2I3.2F6.1.1X.11I4.10H  LOSS..DB)  0057 

59  FORMAT !4X»I3»2F6.1/20X»11< 13. Al) »6H  M0DE/20X , 1 1 1  4 ,7H  ANGLE /4X  .  0058 

2I3.2F6.1  .1X.11I4.5H  DBU )  0059 

24  FORMAT  (6X»4HPWR=F6.2*2HK.W»42X»8HREQ.S/N=I3»2HDB)  0060 

70  FORMAT ( 6X . 3HGMT .2X.3HLUF.3X.3HFOT .3!  5X . 3HGMT . 2X . 3HLUF » 3X . 3HFOT ) / 1  0061 

1H04(4X. I3.6X.F6.1 ) / 1 H04 < 4X . I  3 . 2F6 . 1 ) / 1H04 ( 4X , I  3 ,6X . F6 . 1 ) / 1H04 ! 4X . I  0062 
23.2F6.1)/1H04!4X.I3.6X,F6.1)/1H04(4X,I3.2F6.1)/1H  Al )  0063 

94  FORMAT (8X.I3.2F6.1.1X.11I4)  0064 

907  FORMAT! 1H0/1H0/1H0)  0065 

908  FORMAT ! 1H  )  0066 

909  FORMAT! 1H1 )  0067 

FREL! 11 )=500.  0068 

FREL ! 12 >=500.  0069 

DO  6461  I A= 1 .10  0070 

6461  I FR ( I A ) =FREL ( I  A )  0071 

:  BEARINGS  OF  RHOMB1CS 

IF(XETA)  1320.1320.1321  0072 

1321  XETA=ABSF(XETA-XTR)  0073 

IF(XETA-180. )  1320,1320.211  0074 

211  XETA=360.-XETA  0075 

1320  I F ! WET  A )  1323.1323.1324  0076 

1324  WETA=ABSF(WETA-BTRY)  0077 

I F (WET  A-180 • )  1323.1323.210  0078 
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210  WETA=360.-WETA  0079 

1323  XQ=  *001  0080 

RQ=  #001  0081 

XP=4.  0082 

RP  =  4 •  0083 

I NH=XNH  0084 

I NL=XNL  0085 

I ND=XND  0086 

JNH=RNH  0087 

JNL=RNL  0088 

JND=RND  0089 

N IG=0  0090 

J I G= J I G+l  0091 

MAsX ABSF ( MAN )  0092 

C  CONSTANTS  TO  BE  USED  FREQUENTLY 

TOP(1)*20.  0093 

TOP!  2  )  *20 .  0094 

TOP ( 3 ) =25 •  0095 

TOP ( 4 ) =23. 5  0096 

R0*6376.6  0097 

XIP=1  0098 

I F ( MAP )  5612.5612.5613  0099 

5612  WRITE  OUTPUT  TAPE  6.3.  0100 

1  nocir.bmons.ssn.ha.har.xlat.ax.xlong.al.ylat.ay.ylong.yl.  0101 

2XTR.BTRY.GCDNM  0102 

5613  FK=69 . 057  0103 

LZP= 1  0104 

L I P* 1  0105 

M I P=2  0106 

IF(IHR)  8400.8400.3219  0107 

8400  GO  TO  (3219.3219.3220,3221.3219,3219) .METHOD  0108 

3220  MIP»1  0109 

GO  TO  3219  0110 

3221  MA=0  0111 

3219  DO  3214  I T =M I P , 24 ,M I P  0112 

FREL(ll)-FOTdT)  0113 

JUG=0  0114 

KAT  = 1  0115 

DO  251  IJ=1.11  0116 

NANGLE ( I J ) =0  0117 
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MODE  (  I  J  )  =0  0118 

251  ABC  (  I J  )  =  AZZ  0119 

FREQ  =  FREL  ( 1  )  0120 

DFREO  =  3  •  0121 

NF  =  1  0122 

C  F2  REGION  LAYER  HEIGHTS 

FL  =  0.  0123 

DO  926  1 1  =K1  »K3  0124 

CALL  POLY ( 2  »  7  »7  *  CLCMII.IT)  / 10 , -1 . 2 *CL AT ( 1 1) /90 . * Y )  0125 

926  FL=FL+Y  0126 

RK=K3  0127 

FL  =  FL/RK*100,  0128 

GO  TO  (927,928)  *K2  0129 

C  MODE  CHOOSING  F2  AND  E-LAYER 

928  ME(1 )=GCDKM/2000*+*9  0130 

ME  (  2  )  =ME  ( 1  )  +  l  0131 

MF( 1 )=GCDKM/4000.+.8  0132 

MF  (  2  )  =MF  ( 1  )  +  1  0133 

MF  (  3  )  =MF  (  2  )  + 1  0134 

GO  TO  929  0135 

927  ME  ( 1  )  =1  0136 

ME  (  2  )  =  2  0137 

MF  C 1  )  =1  0138 

MF  (  2  )  =2  0139 

MF  (  3  )  =3  0140 

C  RADIATION  ANGLE  —  E-LAYER  MODES 

929  DO  930  IK=1*2  0141 

I  M=ME  (IK)  0142 

P I  =ME  (  IK)  0143 

ELD(  IM)=GCDKM/PI  '  0144 

EDEL ( I M ) =ATANF ( (COSF(GCD/ ( 2.*PI ) ) -RO/ ( RO+1 10 • ) ) /S INF ( GCD/ ( 2«*PI )) )  0145 

1*BK  0146 

IF ( EDEL ( IM ) )  1032,930,930  0147 

1032  DO  1033  11=1*2  0148 

1033  ME (II) =ME (II )+l  0149 

GO  TO  929  0150 

930  CONTINUE  0151 

C  RADIATION  ANGLE —  F2-LAYER  MODES 

1934  DO  931  I  K  =  1 »  3  0152 

IM=MF(IK)  0153 
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P I  =  I  M  °154 

FLD( IM)=GCDKM/PI  0155 

FDEL ( IM) =ATANF ( ( COSF  <  GCD/ ( 2  .  *P I  ) ) -RO/ ( RO+FL  )  ) /SINF ( GCD/ (2.*PI ) )  )  0156 

1*BK  0151 

I F ( FDE  L  (  I M  )  )  932  ♦931*931  0158 

932  DO  1933  11=1*3  0159 

1933  MF (II) =  MF ( I  I )  +  l  0160 

GO  TO  1934  0161 

931  CONTINUE  0162 

GO  TO<  933*934) *K2  0163 

C  PENETRATION  FREQUENCIES  — E-REGION 

934  DZ=EMF ( 4* I T )  0164 

DY  =  EMF ( 5  » I  T  )  0165 

IF(DZ-DY)  940*940*941  0166 

940  ZR  =  DZ  0167 

12=2  0168 

ZS=DY  0169 

GO  TO  939  0170 

941  ZR=DY  0171 

12=3  0172 

ZS=DZ  0173 

GO  TO  939  0174 

933  ZR  =  EMF(  1  *  I  T  )  0175 

ZS=ZR  0176 

939  DO  1663  11=1.10  0177 

ELL ( II ) =0*  0178 

1663  FLL (  I  I ) =0.  0179 

DO  827  11=1.2  0180 

I m=ME (II)  0181 

ERCR ( I M ) =0  •  0182 

EXMT ( I M ) =0  .  0183 

C  MODE  ELIMINATION  E  AND  F2  REGIONS 

ARZ=ELD( IM)*DK  0184 

I F ( 1 6 • — ARZ )  2900.2900,2901  0185 

2900  ELFK  =  1 • 02  0186 

GO  TO  2902  0187 

2901  ELFK  = ( <  ( < ( <-4.368460907E-9*ARZ  +  1.334494261E-7)*ARZ-5.976618436E-6>  0188 

l*ARZ+2. 62480831 5 E-4 )* ARZ- 5. 0384 76266 E-3 )* ARZ +3. 76 138 5053 E-2 ) *ARZ"1  0189 
2.133200756E-2)*ARZ+.2085  0190 

2902  Y=ZR*ELFK+.05  0191 
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IF(Y-FREQ)  328,328.827 
328  ELL ( IM)=1000. 

827  CONTINUE 

DO  227  11=1,3 
I M=MF (  I  I ) 

FRCR ( I M  )  =0  » 

FXMT ( I M ) =0  , 

ARZ=RO*2 •* ( P I 2-FDEL ( I M ) *AK-AS I NF ( RO*S I NF ( P I 2+FDEL ( I M ) *AK ) / (RO+llO. 
1  )  )  )*DK 

I F ( 16 • -ARZ )  2903,2903,2904 
2903  ELFK= 1 • 02 
GO  TO  2905 

2©04  ELFK= ( ( ( ( ( ( -4, 36 8 46 0907 E ~9* ARZ +1 ,334494261E-7)*ARZ-5,9766l8436E-6) 
l*ARZ+2.624808315E-4)*ARZ-5.038476266E-3 )*ARZ+3.761385053E-2)*ARZ-1 
2ol33200756E~2)*ARZ+,2085 
2905  Y=ELFK*ZS+,05 

IF(Y-FREQ)  227,227,228 
228  FLL ( IM ) =1000 • 

227  CONTINUE 

DO  946  I  I  =  K2  »K3 

F24 ( 1 1 )  =  ( F2S ( I  I » I T )* ( 180,-SSN )+F2H( I  I ♦ I T ) * ( SSN-10 . ) >/170. 

946  GMA ( 1 1  I  =  ! GML ( I  I , I T )*( 180,-SSN )+GMH(  I  I , I T ) * ( SSN- 10 • )  )/170. 

Sl=  (  FREQ-GMA  (  K 2  >  ) / ( F24 ( K2  )  -GMA  (  K. 2  )  ) 

S2=< FREQ-GMA <K3) )  / ( F24 ( K3 ) ~GMA ( K  3 ) ) 

I F ( S1-S2 )  947,947,948 

947  R=S2 

GO  TO  949 

948  R=S1 

949  DO  945  IK=1,3 
I M=MF ( IK) 

I F ( FLL ( IM ) )  508,508,945 
508  ARK=FLD<IM)#DK 

FLF ( I M )=(((( ( (-6,712654756E-9*ARK+4.49151441E-7)*ARK-9,985831104E- 
16 )*ARK+6. 8652598 17E-5)*ARK+9,202437332E-5)*ARK+2.264634341E-3>*ARK 
2+4,699243101E-3)*ARK 
I F ( FLF ( IM ) -R )  950,945,945 

950  FLL ( IM)=1000. 

945  CONTINUE 

GYR=(GY(K2)+GY(K3) )/2. 

GO  TO  ( 1740,9119 ) .KIP 
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C  MAJOR  LAND  BODIES 

1740  WLD=0.  0232 

DO  826  II=K1»K3  0233 

IF  ( CLONG (II))  1702. 1702*1701  0234 

1701  CZG  =  360. -CLONG ( I  I )  0235 

GO  TO  1703  0236 

1702  CZG=ABSF ( CLONG ( II))  0237 

1703  CALL  NOISY ( 6.X2.CZG.Y )  0238 

826  WLD=WLD+Y  0239 

WD=WLD/RK  0240 

9119  IF(WD)  953.954,954  0241 

953  ER=80.  0242 

S I  GMA=  5 •  0243 

GO  TO  955  0244 

954  ER  =  4.  0245 

S  I GMA= .001  0246 

C  E-LAYER  MODES 

C  GROUND  LOSS  REFLECTION  FACTORS  AND  LOSS 

955  DO  726  11=1.2  0247 

I M=ME (II)  0248 

IF  (ELL( IM) )600, 600,726  0249 

600  RS= I M- 1  0250 

CALL  LOSS(EDEL(IM) ,ER. FREQ. SIGMA, Y)  0251 

GE ( I M ) = ABSF ( Y*RS )  0252 

R I = I M  0253 

C  SK.YWAVE  ABSORPTION 

ESK  Y  (  IM)=615.5*RI*ABI  (  I  T  )  /  (  COSF  (  AS  I  NF  (  RO*S  I  NF  (  P  I  2  +  EDEL  ('  I M  )  *AK. )  /  (  RO  0  2  54 

1+100.) ) )*(FREQ+1.12*GYR)**1.98)  0255 

RS= I M  0256 

C  ANGLE  AT  IONOSPHERE 

PHE=RO*SINF(PI2+EDEL( I M ) +AK ) / ( RO+ 110. )  0257 

C  FREE  SPACE  LOSS 

ZOR=  GCD/(2.*RS)  0258 

FR=(2.*RS*SINF(ZOR)*RO/PHE)*EEK  0259 

FSE ( IM)=36.58+20.*LOG10F( FR) +20.*LOG10F( FREQ)  0260 

726  CONTINUE  0261 

C  F2-LAYER  MODES 

C  GROUND  LOSS  REFLECTION  FACTORS  AND  LOSS 

DO  626  11=1.3  0262 

I M=MF (II)  0263 
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I F ( F LL ( I M )  )  60 1 »  60 1  *  626  0264 

601  RS=  I M- 1  0265 

CALL  LOSS(FDEL( IM) ,ER»FREQ»SIGMA,Y)  0266 

GF( IM)=ABSF< Y*RS)  0267 

R I  =  I  M  0266 

C  SK.Y  WAVE  ABSORPTION 

F  SKY ( IM)=615.5*RI*ABI < I  T  )  / ( COSF < AS  I NF < RO*S I NF ( P 1 2+FDEL < IM)*AK)/(RO  0269 

1+100. )) )*<FREQ+1.12*GYR)**1.98)  0270 

RS=  I M  0271 

C  ANGLE  AT  IONOSPHERE 

PHE=RO*S I NF ( P I 2+FDEL ( I M ) *AK ) / ( RO+FL )  0272 

C  FREE  SPACE  LOSS 

ZOR=  GCD/(2.*RS)  0273 

FR= ( 2 • *RS*S I NF ( ZOR ) *RO/PHE )  *EEK  0274 

F  SF ( IM)=36®58  +  20®  *L0G1 OF ( FR ) +20 • *LOG 1  OF ( FREQ )  0275 

626  CONTINUE  0276 

C  ANTENNA  DETERMINATION 

KAP=XA8SF(  IANT  )  0277 

IF(IANT)  18  ,181,181  0278 

181  KAP=4  0279 

DO  82  11*1,2  0280 

I M*ME (II)  0281 

82  EXMT ( I M )  =  I  ANT  0282 

DO  62  11*1,3  0283 

I M=MF (II)  0284 

62  FXMT ( I M )  =  I  ANT  0285 

GO  TO  1185  0286 

18  GO  TO  ( 31,32,33) ,KAP  0287 

31  T  ANT ( 1 ) =ARA ( 1 )  0288 

TANT ( 2 ) =ARA ( 2 )  0289 

GO  TO  180  0290 

32  T ANT ( 1 ) =AVA ( 1 )  0291 

T ANT ( 2  >  =AVA ( 2 )  0292 

GO  TO  180  0293 

33  T ANT ( 1 ) =AHA ( 1 )  0294 

T  ANT ( 2 ) =AHA ( 2 )  0295 

180  DO  182  11*1,2  0296 

I M=ME (II)  0297 

IF(ELL( IM) >603,603,182  0298 

603  CALL  GAIN(ICAP,EOEL(  IM)  ,XETA,FREQ,XQ,XP,XND,XNL,XNH,EXMT(  IM)  )  0299 
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182  CONTINUE  0300 

C  POWER  GAIN  OF  TRANSMITTING  ANTENNA 

DO  382  11=1*3  0301 

I M=MF (II)  0302 

I F ( FLL ( I M ) )  605*605,382  0303 

605  CALL  GAIN( KAP,FDEL ( IM) ,XETA»FREQ»XQ»XP ,XND»XNL .XNh.FXMT ( IM) )  0304 

382  CONTINUE  0305 

1185  I F (METHOD-4 )  185,184*185  0306 

184  I F ( I  HR )  185*  185*  183  0307 

183  I ANR  =0  0308 

MA=0  0309 

185  KUP=XABSF( IANR )  0310 

IF  (IANR)  187*281*281  0311 

281  KUP  =  4  0312 

DO  482  11=1.3  0313 

I M=MF (II)  0314 

482  FRCRt  I M  )  =  I ANR  0315 

DO  582  11=1*2  0316 

I M=ME (II)  0317 

582  ERCR ( I M ) = I ANR  0318 

GO  TO  1285  0319 

187  GO  TO  (741*42*43) *KUP  0320 

741  RANT ( 1 ) =AR A ( 1 )  0321 

RANT ( 2 ) =ARA ( 2 )  0322 

GO  TO  280  0323 

42  RANT ( 1  )  *AVA ( 1)  0324 

RANT ( 2  )  =AVA ( 2 )  0325 

GO  TO  280  0326 

43  RANT ( 1 )=AHA( 1 )  0327 

R ANT ( 2 ) = AHA ( 2 )  0328 

C  RESPONSE  OF  RECEIVING  ANTENNA 

280  DO  282  11=1.3  0329 

I M=MF ( I  I )  0330 

I F ( FLL ( I M ) )  607*607,282  0331 

607  CALL  GAIN(KUP*FDEL( IM) ,WETA*FREQ»R0»RP»RND*RNL*RNH*FRCR( IM) )  0332 

282  CONTINUE  0333 

DO  782  11=1.2  0334 

IM=ME(II)  0335 

I F ( ELL ( IM) )609*609*782  0336 

609  CALL  GAIN(KUP,EDEL( IM) ,WETA»FREQ.RQ*RP,RND.RNL,RNh,ERCR( IM) )  0337 
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782  CONTINUE  0338 

1285  I F ( MAP )  285 ♦ 285  *  388  0339 

285  GO  TO  ( 286.388) .LIP  0340 

286  GO  TO  ( 330.330.330,1331 ) »KAP  0341 

1331  GO  TO  (232.232.232.233) »KUP  0342 

233  WRITE  OUTPUT  TAPE  6 , 1 3 . I  ANT , MA , I ANR  0343 

GO  TO  388  0344 

330  GO  TO  (335.335,335,336) ,KUP  0345 

335  WRITE  OUTPUT  TAPE  6.1,  0346 

1  TANT.INH.INL.IND.MA.RANT, JNH, JNL.JND  0347 

GO  TO  388  0348 

336  WRITE  OUTPUT  TAPE  6.2,  0349 

1  TANT.INH.INL.JND.MA.IANR  0350 

GO  TO  388  0351 

232  WRITE  OUTPUT  TAPE  6 .4 , I  ANT  ,MA .R ANT ♦ JNH . JNL . JND  0352 

388  L I P=2  0353 

DO  333  II=K2»K3  0354 

T I P=0 •  0355 

C  DETERMINATION  OF  TYPE  OF  CIRCUITtPOLAR, TEMPERATE, ETC) 

RM=ABSF(GLAT( I  I ) )  0356 

IF(RM-60.*AK)  331,331,332  0357 

331  I D= 1  0358 

IF(ABIY( 1,  IT  )  )  9732.9732.333  0359 

9732  T  I  P  =  5 •  0360 

GO  TO  333  0361 

332  IF (RM-70.*AK  )  334,334,  35  0362 

334  I D=3  0363 

IF(K2-1>  333,333,  36  0364 

36  I D=4  0365 

GO  TO  1725  0366 

35  I D=2  0367 

333  CONTINUE  0368 

C  TRANSMISSION  LOSS  E-LAYER  MODES 

1725  DO  960  11*1,2  0369 

I M=ME (II)  0370 

960  ELL  (  IM)=ESKY(IM)  4-GE  (  I M  )  +FSE  (  I M )  -EXMT  (  IM)-ERCR(  IM)+ADJ ( ID)+ELL( IM)  0371 

LD=ME ( 1 )  0372 

LE*ME ( 2 )  0373 

I F ( ELL ( LD ) -ELL ( LE ) )  961.961,962  0374 

961  ELOS=ELL ( LD )  0375 
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MDE=LD  0376 

GO  TO  963  0377 

962  ELOS=ELL(LE>  0378 

MDE=LE  0379 

C  TRANSMISSION  LOSS  F2-LAVER  MODES 

963  DO  964  11=1.3  0380 

I M=MF (II)  0381 

964  FLU  I M ) =FSKY ( IM)+GF( IM)+FSF( IM)-FXMT(  IM)-FRCR( IM)+ADJ( ID)+FLL( IM)  0382 

LG=MF ( 1 )  0383 

LH=MF ( 2 )  0384 

L J=MF ( 3 )  0385 

IF(FLL(LG)-FLL(LH) )  965.965.966  0386 

965  FL1=FLL ( LG )  0387 

MDF=LG  0388 

GO  TO  967  0389 

966  FL1=FLL ( LH )  0390 

MDF=LH  0391 

967  I F ( FL1-FLL ( L J ) )  968,968.969  0392 

968  FL0S=FL1  0393 

GO  TO  970  0394 

969  FLOS=FLL(LJ)  0395 

MDF=L J  0396 

970  I F ( ELOS-FLOS )  971,971,972  0397 

971  NXLOS ( NF ) =ELOS  0398 

N ANGLE (NF)=EDEL(MDE)+«5  0399 

MODE ( NF ) =MDE  0400 

ABC ( NF ) =AEC  0401 

GO  TO  911  0402 

972  NXLOS(NF)=FLOS  0403 

N ANGLE ( NF ) =FDEL ( MDF ) +, 5  0404 

MODE ( NF ) =MDF  0405 

ABC ( NF ) =AFC  0406 

C  MIXED  MODES  OF  PROPAGATION 

911  FHOP*MF ( 1 ) - 1  0407 

IF(FHOP)  3903.3903,3898  0408 

3898  XHOP=FHOP+1.  0409 

C  AVERAGE  HEIGHT  OF  REFLECTING  LAYERS 

EFHT  =  (  HO.+FHOP*FL  )  /  (  FHOP  +  1.  )  0410 

C  RADIATION  ANGLE 

EFDEL=ATANF! (COSF(GCD/(2.*XHOP) )-RO/ ( RO+EFHT ) ) /SI NF ( GCD/ ( 2 .*XHOP ) ) 


0411 
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1)*BK  0412 

C  ANGLE  AT  IONOSPHERE 

PHE=RO*SINF(PI2+EFDEL#AK) /(RO+llO. )  0413 

ARZ=79. 138*( PI 2-EFDEL*AK-ASINF ( PHE ) )  0414 

I F ( EFDEL-1  •  )  3899*3900*3901  0415 

3899  FHOP=FHOP+l.  0416 

GO  TO  3898  0417 

C  MODE  ELIMINATION 

3900  ELFK=1.02  0418 

GO  TO  3902  0419 

3901  ELFK=( <  (  I  I <-4.368460907E-9*ARZ+1.334494261E-7)*ARZ-5.976618436E-6>  0420 

1*ARZ +2. 62480831 5 E-4 ) *ARZ- 5.0 38476266 E-3 )*ARZ +3. 76 138 505 3E-2 ) *ARZ~1  0421 
2.133200756E-2>*ARZ+.2085  0422 

3902  Y  =  ZS*E  LFK  0423 

4332  IF(Y-FREO)  3903.3905*3905  0424 

3903  NOSINF ) =1000.  0425 

GO  TO  612  0426 

3905  Y  =ZR*ELFK  0427 

IF(Y-FREQ)  650.650*651  0428 

651  NOSINF )=1000.  0429 

GO  TO  612  0430 

650  ARK= (GCDKM*DK-ARZ ) /FHOP  0431 

FLFK  =(<(((  I-6.712654756E-9*ARK  +  4.49151441E-7)*ARK-9.985831104E-  0432 

16  )*ARK  +  6. 8652 5981  7E-5  )  *ARK+9 . 20 24 3 73 3 2E -5  )  *ARK+2.26463434lE-3  )*ARK  0433 

2+4.699243101E-3 ) *ARK  0434 

REF= ( FREQ-GMA (I2))/(F24(I2  J-GMA (12))  0435 

IF  (FLFK.-REF  )  652.653.653  0436 

652  NOSINF )=1000.  0437 

GO  TO  612  0438 

C  SKY  WAVE  ABSORPTION 

653  EFSKY  =  615.5*XHOP*ABIY ( 1 , I T ) / ( COSF ( AS  I NF ( RO*S INF (PI2  +  EFDEL*AK ) / ( RO  0439 

1  +  100. )  ) ) * ( FREQ+1 . 12*GYR  >**1.98 )  0440 

C  GROUND  REFLECTION  LOSS 

3908  CALL  LOSS ( EFDEL . ER »FREQ ♦ S I GMA , Y )  0441 

3909  GEF=ABSF (Y*FHOP)  0442 

C  ANGLE  AT  IONOSPHERE 

SPHE=RO*SINF ( P I 2+EFDEL+AK ) / ( RO+EFHT )  0443 

C  FREE  SPACE  LOSS 

ZOR=GCD/ ( 2 •* I FHOP+ 1 . ) )  0444 

FR=(2.*«FHOP+l. )*SINF(ZOR)*RO/SPHE)*EEK  0445 
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EFSF=36. 58+20.*L0Gl0F ( FR ) +20 . *LOG 10F ( FREQ )  0446 

GO  TO  ( 3180,3180,3180,3181 ) »KAP  0447 

3181  EFXMT  =  I  ANT  0448 

GO  TO  3182  0449 

C  POWER  GAIN  OF  TRANSMITTING  ANTENNA 

3180  CALL  GAINtKAP, EFDEL,XETA,FREO,XQ,XP,XND,XNL,XNH, EFXMT )  0450 

3182  GO  TO  (3183, 3183, 3183, 3184), KUP  0451 

3184  EFRCR=IANR  0452 

GO  TO  3190  0453 

C  RESPONSE  OF  RECEIVING  ANTENNA 

3183  CALL  GAIN( KUP ,EFDEL,WETA,FREO,RQ,RP ,RND,RNL ,RNH ,EFRCR )  0454 

C  TRANSMISSION  LOSS 

3190  NOS(NF)=EFSKY+GEF+EFSF-EFXMT-EFKCR+ADJ( ID)  0455 

612  IF(NOS(NF)-NXLOS(NF) )  623,622,622  0456 

623  NANGLE(NF)=EFDEL+,5  0457 

MODE (NF)=FHOP+l,  0458 

B  ABC(NF  1=676060606060  0459 

NXLOSI NF) =NOS( NF )  0460 

622  IF(NXLOS(NF)-990  )  6911,624,624  0461 

624  ABC ( NF ) = A2 Z  0462 

NANGLE ( NF ) =0  0463 

MODE ( NF ) =0  0464 

GO  TO  8740  0465 

6911  KA J= 1  0466 

GO  TO  ( 6912,6912,6912,9779,6912  ,6912 )  .METHOD  0467 

6912  POR=10.*LOG10F(PWR*10QO. )  0468 

I  F ( JUG) 3869,3869,8712  0469 

C  1MC/S  ATMOSPHERIC  NOISE  DETERMINATION 

3869  CC=CKC ( IT )  0470 

C  SEASON  AND  HOUR  BLOCK  FOR  ATMOSPHERIC  NOISE  DETERMINATION  0471 

L I B=4  0472 

I  CC  =  CC/2 •  0473 

I  F ( ICC)  8889,8889,8888  0474 

8888  GO  TO  (8701,8701,8702.8702.8705,8705,8706,8706,8708,8708,8887,  0475 

18887), ICC  0476 

8887  CC=CC-24.  0477 

8889  K J= 1  0478 

TM=-2 •  0479 

JK=1  0480 

GO  TO  8703  0481 
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8701  K J  =  1  0482 

TM=2 •  0483 

JK  =  5  0484 

GO  TO  8703  .  0485 

8702  K J=5  0486 

TM=6 •  0487 

JK.  =  4  0488 

KAJ=2  0489 

GO  TO  8703  0490 

8705  KJ  =  4  0491 

TM=10.  0492 

JK  =  3  0493 

L 16=5  0494 

GO  TO  8703  0495 

8706  K  J  =  3  0496 

TM=14.  0497 

JK=2  0498 

L 1 6=5  0499 

GO  TO  8703  0500 

8708  KJ  =  2  0501 

TM= 18  •  0502 

JK= 1  0000 

<AJ=2  0503 

8703  IF(Y2>  9702.9702.9701  0504 

9701  CEG=360.-Y2  0505 

GO  TO  9703  0506 

9702  CEG=A8SF(Y2)  0507 

9703  CALL  NOISY(KJ.X2  .CE6.ATN0)  0508 

CALL  NOISY! JK.X2.CEG.ATNY)  0509 

ATNQ=(ATNO+< ATNY-ATNO ) * ( CC-TM ) /4. )  0510 

L06=0  0511 

JUG= JUG+1  0512- 

GO  TO( 8712 .8713 ) .KAJ  0513 

C  FREQUENCY  DEPENDENCE  OF  ATMOSPHERIC  NOISE 

8713  MOT  =MOUSE/ 3  0514 

IF(MOT)  8714.8715.8714  0515 

8715  MOT  =  4  0516 

8714  I F ( X  2 )  8716.8717.8717  0517 

8716  L08= 1  0518 

8717  I F ( MOT-3 )  8718.8719.8719  0519 
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8718  L I B=  5-LOB  0520 

GO  TO  8712  0521 

8719  L  I B=4+L0B  0522 

8712  IF(LIB~4)  2179.2179,2178  0523 

2179  IF(ATNQ-20.)  2130,2178.2178  0524 

2130  ATNQ=20.  0525 

2178  CALL  POLY(LIB,7,7,ATNQ/lO.-5.,FREQ/10.-1.35,y )  0526 

ATNO=Y*10.+130.  0527 

C  GALACTIC  NOISE 

I F ( FREQ-GMA ( 2 ) )  8720,8720,8721  0528 

8721  GNOS  =  ((((((((  ,l095032130E-10*FREQ+.2442852795E-8 )*FREQ-. 15854096  0529 

108E-5)*FREQ+.1513740543E-3)*FREQ-.6306642189E-2 ) *FREQ+. 1390178355E  0530 

2+0)*FREQ-.170l088795E+l )*FREQ+. 1 1 2 1 1 3 0396E+2 ) *FREQ~. 3426395658 E+2 )  0531 

3*FREQ+.2020285454E+3  0532 

GO  TO  9797  0533 

8720  GNOS  =  1000.  0534 

C  MAN  MADE  NOISE 

9797  I F ( MAN )  872 2 , 872 3 , 872 3  0535 

8722  KJ=XABSF(MAN)+5  0536 

S0B=0.  0537 

CALL  POLY  (KJ, 1,10, SOB, FREQ. Y)  0538 

XNO I S= Y  0539 

GO  TO  8724  ff540 

8723  XNO I S=MAN  0541 

C  DETERMINATION  OF  CONTROLLING  NOISE 

8724  RCNSE= ATNO  0542 

I F { ATNO-GNOS)  8725,8725,8726  0543 

8726  RCNSE=GNOS  0544 

8725  IF(RCNSE-XNOIS)  8729,8729,8728  0545 

8728  RCNSE=XNOIS  0546 

C  AVAILABLE  S I GNAL-TO-NO I SE  RATIO 

8729  XLOS=NXLOS(NF  )  0547 

ROT  =  XLOS-POR  0548 

GOT  =  RCNSE-ROT  0549 

GO  TO( 505,505.505,505,505.506) .METHOD  0550 

506  NXLOS ( NF ) =GOT  0551 

GO  TO  9777  0552 

505  WANT  =GOT-RSN  0553 

C  CIRCUIT  RELIABILITY 

IF(ABSF(WANT)-(TOP( I D ) +T I P ) ) 867 ,867 . 868 


0554 
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868  REL=5IGNF(100.»WANT)  0555 

GO  TO  869  0556 

867  K.J=  ID+10  0557 

GO  TO  <888. 8733*8733. 8733)  »IO  0558 

888  IF(ABIY< 1.  IT  )  )  8733.8733.8732  0559 

8732  KJ=KJ-1  0560 

8733  CALL  POLY < K J » 5 • 1 0 . FREQ/ 10.-1 . 35 . WANT/ 1 0 •  . Y )  0561 

REL  =  Y*  1 0  e-«-50  »  0562 

869  IF  I  ID— 1 )  8734.8734.8735  0563 

8735  XRL=90.  0564 

GO  TO  8736  0565 

8734  XRL=99.  0566 

8736  IF(REL-XRL)  8739.8739.8738  0567 

8738  REL=XRL  0568 

8739  IF(REL-1.  )  8740.8740,8741  0569 

8740  REL=0.  0570 

8741  NXLOS ( NF ) *REL  0571 

K I  P  =  2  0572 

I F ( METHOD-3 )  9777 , 3207 e 9777  0573 

C  ITERATION  FOR  LUF 

3207  I F ( REL-90. J  3201.3202.3203  0574 

3201  GO  TO  ( 3204.3202 )  ,KAT  0575 

3202  FF ( I T ) =FREQ  0576 

GO  TO  3214  0577 

3203  I F ( FREQ-3. )  3208.3208,3209  0578 

3208  FF ( I T ) = — 3 •  0579 

GO  TO  3214  0580 

3209  FREQ=FREQ-.1*DFREQ  0581 

IF(FREO-3. >3208.4850,4850  0582 

4850  K AT=2  0583 

GO  TO  939  0584 

3204  Z J=K J  0585 

GO  TO  9777  0586 

9779  I F ( I  HR )  9777,9777.9778  0587 

C  FIELD  STRENGTH 

9778  A LOSS =— NXLOS ( NF )  0588 

NXLOS ( NF ) =ALOSS+107.2  +  8.6859*LOGF ( FREQ ) +4. 343*LOGF  < 1000. *PWR )  0589 

C  FURTHER  ITERATION  FOR  LUF 

9777  NF=NF+1  0590 

IF(NF-ll)  914,914,260  0591 
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914  FREQ=FREL(NF)  0592 

I F ( FREQ-FREL (11))  917,917,916  0593 

917  I F ( FREQ-30, )  939,939,916  0594 

916  NXLOS ( NF ) =0  0595 

FF ( I T ) =0«  0596 

GO  TO  9777  0597 

260  GO  TO  (250,250,3214,3212*7007,250) .METHOD  0598 

7007  GO  TO  ( 7009,7008 ) ,LZP  0599 

7009  WRITE  OUTPUT  TAPE  6 , 34  ,PWR  ,  I RSN  ,  I FR  0600 

7008  WRITE  OUTPUT  TAPE  6,94,  0601 

1  IG< IT) ,UF( IT ) , FOT ( IT), NXLOS  0602 

NIG=NI G+l  0603 

GO  TO  ( 3214,3214,21 ), NIG  0604 

21  WRITE  OUTPUT  TAPE  6,908  0605 

N IG=  0  0606 

GO  TO  3214  0607 

250  GO  TO  ( 2250, 2666 ) ,LZP  0608 

2250  WRITE  OUTPUT  TAPE  6 , 14  ,PWR  ,  I RSN ♦ I FR  0609 

2666  GO  TO  ( 2667 , 2667 , 2667 , 266 7 , 2667 ♦ 2668 ), METHOD  0610 

2668  WRITE  OUTPUT  TAPE  6,29,  0611 

1  IG( IT-1 ) ,UF( IT-1 ) ,FOT ( IT-1 ) , ( MODE (  I  )  ,A0C(  I ) .1=1 ,11 ) .NANGLE  0612 

2  * IG (  IT)  *UF  <  IT)  ,FOT( IT)  ,  NXLOS  0613 

GO  TO  3214  0614 

2667  WRITE  OUTPUT  TAPE  6,9,  0615 

1  IG( IT-1) ,UF( IT— 1 ) »FOT ( IT— 1) *( MODE ( I ) »ABC( I ) » I =1 , 11 ) .NANGLE  0616 

2»IG( IT)  ,UF( IT) ,FOT( IT) , NXLOS  0617 

GO  TO  3214  0618 

3212  GO  TO  ( 3213.3215 ) ,LZP  0619 

3213  IF(IHR)  8002,8002.8003  0620 

8003  WRITE  OUTPUT  TAPE  6,80,PWR.IFR  0621 

GO  TO  77  0622 

8002  WRITE  OUTPUT  TAPE  6,8, IFR  0623 

3215  IF(IHR)  69,69,77  0624 

69  WRITE  OUTPUT  TAPE  6,19,  0625 

1  IG ( IT-1) ,UF(  IT-1)  , FOT ( IT- 1) ,<MODE<  I ) ,ABC( I ) , I = 1 » 1 1 ). NANGLE  0626 

2  » IG( IT) *UF (IT) ,FOT ( IT) .NXLOS  0627 

GO  TO  3214  0628 

77  WRITE  OUTPUT  TAPE  6,59,  0629 

1  IG( IT-1) *UF< IT-1 > ,FOT( IT-1) »<MODE( I ) ,ABC< I ) . I = 1 ♦ 1 1 ) . NANGLE  0630 

2  ,  IG ( IT)  »UF < I T )  »FOT ( IT ) .NXLOS  0631 
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3214  LZP  =  2  0632 

GO  TO  (2270. 51*290.51. 41. 51) .METHOD  0633 

41  GO  TO  (40.51). JIG  0634 

40  WRITE  OUTPUT  TAPE  6.907  0635 

GO  TO  270  0636 

51  WRITE  OUTPUT  TAPE  6.909  0637 

J I G=0  0638 

GO  TO  270  0639 

290  I F ( MAP  )  6228,6228.6229  0640 

6229  CALL  CURVY ( FOT  .FF )  0641 

GO  TO  270  0642 

6228  WRITE  OUTPUT  TAPE  6 .24 ,PWR . I RSN  0643 

IF(IHR)  2291.2291.2290  0644 

2291  WRITE  OUTPUT  TAPE  6.7,  0645 

1  ( IG( I ) ,FF( I ) ,FOT( I ) ,IG( 1+6 ) .FF( 1+6) .FOT ( 1+6) »IG( 1  +  12 )  ,  0646 

2  FF ( 1+12) .FOT ( 1+12) . IG( 1+18) .FF( 1+18) »FOT( 1+18 ) .1=1.6) »AAA  0647 

GO  TO  2270  0648 

2290  WRITE  OUTPUT  TAPE  6,70,  0649 

1  ( I G ( I ) ,FOT( I ) , I G ( 1+6 ) .FOT( 1+6) , IG( 1+12) ,FOT( 1+12 ) . IG( 1+1  0650 

*  \  ' 

28 ) .FOT ( 1  +  18) , I G ( 1+1 ) »FF( 1+1 ) ♦  FOT( 1  +  1  )  ,IG( 1  +  7)  ,FF( 1+7 ) ,FOT ( I +7 ) » IG  0651 

3 ( 1+13 ) *FF( 1  +  13 ) »F0T ( I +1 3 >  » I G < I +  19 ) ♦ FF ( I +  19 ) . FOT ( I +  19 ) . I  =  1 .6 ♦ 2 ) » AAA  0652 

2270  I F ( 3-MI T )  1969.1969.1970  0653 

1969  WRITE  OUTPUT  TAPE  6,909  0654 

MI  T  =  1  0655 

GO  TO  270  0656 

1970  MIT=MIT+1  0657 

270  RETURN  0658 

END  0659 
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SUBROUTINE  CURVY (COT* BOT)  0001 

GENERATES  LINE  GRAPHS 

DIMENSION  C0T(24) »ROT (73) »BOT(24) *WOT(73)  0002 

DIMENSION  A0T(73) ,COB( 55) ,BOB(55) *XOB(73)  0003 

DIMENSION  FDEK(IO) »EDEK(10)  0004 

DIMENSION  A MON ( 12 ) *GLAT( 5) »ABI (24)»ABIY(5*24) *CLCK( 5*24) »GY( 5) ♦  0005 

1F2S( 5*24) »F2H( 5*24) ,GML( 5*24) »GMH( 5*24) *MF(3) »ME( 2) *ELL (10)*  FLL ( 10  0006 

2 ) »ELD ( 10 ) , EDEL ( 10 ) *FLD( 10) *FDEL ( 1 0 ) »F24( 5) *GMA (5) * FLF ( 1 0 ) ♦ FSKY ( 10 )  0007 

3  »ESKY( 10) ,GE( 10) * GF ( 1 0 ) , FSE ( 1 0 ) *FSF( 10) »NOS( 11 ) *FXMT( 10 ) *EXMT( 10) *  0008 

4FRCR (10) ,ERCR(10)  » AD J ( 4 ) * ABC ( 1 2 ) *  MODE ( 12 )  0009 

DIMENSION  Q( 20*60*4) *IL(4) » JL ( 4 ) »KL(4) *LK ( 4 ) »JAL(4)  *CKC(24)  0010 

1*G(60)*AB( 60) *5(20*24) ♦ C ( 20 » 24 ) » GAMMA ( 5  *  24 ) »BA ( 60 )  0011 

DIMENSION  RD(5)*SUN(12) » A ( 10  *  7 » 14 ) ♦ RASSN ( 12 ) * EMF ( 5 » 24 )  0012 

DIMENSION  CLAT(5) »CLONG(5) »EMUFY(5»24)*FMUFY(5*24).  UFY ( 5 ) * FOTY ( 5 )  0013 

5  » I G ( 24 )  0014 

DIMENSION  UF(24) »F0T{24) »FF(24) *P(29*16*6) »ABP(2»6)  0015 

DIMENSION  AHA( 2) »ARA(2 ) »AVA(2 ) »RANT( 2 ) *TANT( 2)  0016 

D I  MENS  I ONNXLOS (11) »NANGLE (11)  0017 

DIMENSION  T0P(4) *FREL( 12)  0018 

COMMON  ABP*P  0019 

COMMON  A,Q*G*AB*S,C* GAMMA *BA ♦ EMUFY » FMUFY * IG *UF * FOT *SUN * AVA * ARA *GLA  0020 

IT  * ABI »  ABI Y  *CLCK*  GY  »F2S  *  F2H*GML  *GMH  »MF  *ME  *ELD  *  FLD*  F24, GMA  * ADJ , ABC •  0021 

2 EMF  »CKC»CLAT  *CLONG* AK  *BK*CK*DK»PI2*EK*EEK*GLT  *GLG  » AAA  * ALA  *  AS A ,  0022 

3AFC*ANC,AWC»AEC.ASC*AZZ»MAN,RSN,KW»PWR, I  ANT* I ANR , Y2 *X2 *  MOUSE *SSN *  0023 

4N0CIR»K4*  AX*XL0NG*AL*YLAT*XLAT.YL0NG»YL*BTRY»GCDKM,GCDNM*K1,K2»X  0024 

53»K5*GCD* IRSN»MON,ID*MOUSE*AY»METHOD*AHA*BMONS*MAP  0025 

COMMON  FREL  0026 

COMMON  XNH.XNL*XND,RNH*RNL  »RND » HA  * J I G *HAR ♦ XTR . XETA *WETA * IHR *M IT  0027 

COMMON  KAP»KUP*MA,TANT*INH»INL» IND ♦ RANT  * JNH , JNL » JND  0028 

73  FORMAT? 10X,4HANT= I3*2HDB»20X» 6HNOISE=I 3* 20X*4HANT=I3*2HDB)  0029 

91  FORMAT (6X»2A6»1X*I3,1HH»1X*I3*1HL»1X»I3  »3HDEG  *2X»6HNOISE=I3»2A6.1X  0030 

1 »I3*1HH,1X»I3*1HL*1X*I3»3HDEG)  0031 

72  FORMAT (6X*2A6> IX, I  3 » 1HH ♦ 1 X ♦ I  3 » 1HL » IX , 1 3  * 3HDEG » 1 2X * 6HNOI SE= 1 3  * 12X *  0032 

14HANT=I3»2HDB)  0033 

74  FORMAT ( 10X ,4HANT= 1 3 ,2HDB * 13X . 6HNO ISE= 1 3 » 7X * 2A6 * IX » 13 » 1HH *  IX *  1 3 » 1HL  0034 

1 *1X* I3»3HDEG)  0035 

24  FORMAT ( 10X  »4HPWR  =  F6 • 2*2HKW*42X.8HREQ*S/N  =  I3  *  2HDB )  0036 

3  FORMAT ( 14X* I 5»10X,A6,  10X , 4HSSN  =  F5 , 0  * 10X * A2 » F7, 3 / 1 3X .  11HTR  0037 

2ANSMITTER,13X,8HRECEIVER*12X,8HBEARINGS*6X*7HN.MILES/11X,  F6.2*A1 *  0038 

32H  -  ,F7.2,A1,5X,F6.2,A1»2H  -  * F7 . 2  *  A 1  * 3X * 2F6 . 1 »4X , F8 *  1  )  0039 
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2ANSMI TTER. 13X.8HRECE I VER  .  12X  » 8HBEAR I NGS »6X , 7HN.MI LES/ 1 IX ,  F6.2  ,A1 ,  0038 

32H  -  ,F7.2 »A1 .5X.F6.2 ,A1 »2H  -  * F 7 • 2 ♦ A1 * 3X » 2F6. 1 ,4X » F8  •  1  )  0039 

15  FORMAT! 7X.73A1/6X.74H00  02  04  06  08  10  12  14  0040 

1  16  18  20  22  00/34X . 19HGREENWICH  MEAN  TIME)  0041 

2  FORMAT ( 2X.A1 , IX. A2 .A1 . 73A 1 , A1 , A2 ,  IX  •  A 1 )  0042 

25  FORMAT! 1H1)  0043 

1  FORMAT (1H  )  0044 

WRITE  OUTPUT  TAPE  6.3.  0045 

1  NOCIR.BMONS.SSN.HA.HAR.XLAT .AX.XLONG.AL.YLAT .AY.YLONG.YL.  0046 

2XTR.BTRY.GCDNM  0047 

I  F (METHOD-3 )  4044.4045.4045  0048 


4045 

GO  TO 

( 330,330.330,1331 ) .KAP 

0049 

1331 

GO  TO 

(  232.232.232.233 ) .KUP 

0050 

233 

WRITE 

output  tape  6,73. iant.ma. ianr 

0051 

GO  TO 

388 

0052 

330 

GO  TO 

( 335,335,335,336 ) .KUP 

0053 

335 

WRITE 

OUTPUT  TAPE  6,91, 

0054 

1 

TANT , INH , INL , I ND » MA • KAN T , JNH , JNL , JND 

0055 

GO  TO 

388 

0056 

336 

WRITE 

OUTPUT  TAPE  6,72. 

0057 

1 

TANT, 1NH.INL.IND.MA, IANR 

0058 

GO  TO 

388 

0059 

232 

WRITE 

OUTPUT  TAP£  6.74, 

0060 

1 

I  ANT .MA.RANT , JNH , JNL  » JND 

0061 

388 

L  I  P  =  2 

0062 

WRITE 

OUTPUT  TAPE  6 ,24  ,PWR  ,  I RSN 

0063 

4044 

J  =  0 

0064 

DO  303 

1=1.24 

0065 

X=COT ( 1-1 ) 

0066 

IF(I-l)  202.203.202 

0067 

203 

X=COT(  24) 

0068 

ROT ( 1  ) 

=  X 

0069 

202 

XX  =  COT <  I  ) 

0070 

DO  303 

11=1,3 

0071 

J=J+1 

0072 

CEG=  I  I 

0073 

303 

ROT  (  J+l  ) =X+( XX-X )*CEG/3. 

0074 

I F (METHOD-3 )  6402.6405.6402 

0075 

6405 

J  =  0 

0076 

DO  808 

1=1  .24 

0077 
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X  =  XX 

0078 

IF(I-l)  802.803.802 

0079 

803 

X=B0T ( 24) 

0080 

I  F ( X  )  51.53.54 

0081 

51 

X  =  1  * 

0082 

GO  TO  54 

0083 

53 

X  =  COT<  24) 

0084 

54 

WOK  1  )  =X 

0085 

802 

XX  =  BOT (  I  ) 

0086 

I F ( XX 1  55.57,58 

0087 

55 

XX  =  1  . 

0088 

GO  TO  58 

0089 

57 

XX  =  COT (  I  ) 

0090 

58 

DO  808  11=1,3 

0091 

J  =  J+1 

0092 

CEG= I I 

0093 

808 

WOT ( J+l ) =  X+<  XX-X ) *CEG/3. 

0094 

6402 

DO  608  1=2.72.1 

0095 

B  608 

XOB ( I ) =406060606060 

0096 

DO  20  1=1.73.3 

0097 

B  20 

XOB ( I >=206060606060 

0098 

DO  30  1=1.55 

0099 

B 

BOB ( I ) =606060606060 

0100 

B  30 

COB ( I ) =606060606060 

0101 

B 

COB ( 1 ) =  600360606060 

0102 

B 

COB ( 5 ) =  600460606060 

0103 

B 

COB ( 9 ) =  600560606060 

0104 

B 

COB ( 13 ) =600660606060 

0105 

B 

COB ( 17 ) =600760606060 

0106 

B 

COB ( 21 ) =601060606060 

0107 

B 

COB ( 23 ) =601160606060 

0108 

B 

COB ( 25 1=011260606060 

0109 

B 

COB<  27 ) =010160606060 

0110 

B 

COB ( 29 ) =010260606060 

0111 

B 

COB( 31 1=010360606060 

0112 

B 

COB( 33 ) =010460606060 

0113 

B 

COB ( 35 ) =010560606060 

0114 

B 

COB( 37 ) =010660606060 

0115 

B 

COB ( 391=010760606060 

0116 

B 

COB ( 41 1 =011060606060 

0117 
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B 

COB( 43 ) =011160606060 

0118 

B 

COB<  45  )  =021260606060 

0119 

B 

COB (47 1=020260606060 

0120 

B 

COB( 49 ) =020460606060 

0121 

B 

COB ( 51 ) =020660606060 

0122 

B 

COB( 53 ) =021060606060 

0123 

B 

COB ( 55 ) =031260606060 

0124 

B 

BOB ( 29 ) =456060606060 

0125 

B 

BOB( 21 ) =622060606060 

0126 

B 

BOB( 22 ) =616060606060 

0127 

B 

BOB ( 23 1=236060606060 

0128 

B 

BOB ( 24 1 =446060606060 

0129 

B 

BOB ( 27 1=706060606060 

0130 

B 

BOB ( 28 1 =236060606060 

0131 

B 

BOB ( 30 1 =256060606060 

0132 

B 

BOB ( 31 1 =646060606060 

0133 

B 

BOB (32  1=506060606060 

0134 

B 

BOB ( 33 1 =256060606060 

0135 

B 

BOB ( 34 1=516060606060 

0136 

B 

BOB ( 35 1 =266060606060 

0137 

B 

DOT =546060606060 

0138 

ZOB=3 1 • 

0139 

B 

008=406060606060 

0140 

DO  200  I K= 1 »  55 

0141 

IM=56-I K 

0142 

DO  100  1=1.73 

0143 

B  100 

AOT ( I 1 =606060606060 

0144 

IF(ZOB-2Q.)  4.4,63 

0145 

63 

WOB= 1 • 

0146 

GO  TO  5 

0147 

4 

IF(ZOB-8.)  4062.4062.4063 

0148 

4063 

WQB= • 5 

0149 

GO  TO  5 

0150 

4062 

WOB=.25 

0151 

5 

ZOB=ZOB-WOB 

0152 

I  T  =  1 

0153 

12 

I  SPOT  =  I  T 

0154 

I F ( ROT ( ITl-ZOB)  6.8,8 

0155 

8 

SOB=ZOB+WOB 

0156 

I F ( ROT (IT) -SOB  1  10,6.6 

0157 
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10  AOT (  I  SPOT ) =DOT  015* 

I  T= I T  +  2  0159 

GO  TO  16  0160 

6  I T= I T+ 1  0161 

16  I  F ( I  SPOT-73 )  12.13.13  0162 

13  IF (METHOD-3 )  200. 6565. 20U  0163 

6565  I T= 1  0164 

33  I  SPOT  = I T  0165 

I F ( WOT ( I T ) -200 )  66.38.38  0166 

38  SOB=ZOB+WOB  0167 

I  F ( WOT (IT) -SOB )  80,66.66  0168 

80  IF( AOT( ISPOT l-DOT )  81,66,81  0169 

B  81  AOT ( ISPOT ) =336060606060  0170 

I  T= I T  +  2  0171 

GO  TO  67  0172 

66  I T= I T+l  0173 

67  I F ( I  SPOT-73 )  33.200,200  0174 

200  WRITE  OUTPUT  TAPE  6,2.  0175 

1  BOB(  IM)  .COB(  IM)  .DOB.AOT, DOB. COBUM)  »BOB(  IM>  0176 

14  WRITE  OUTPUT  TAPE  6.15.XOB  0177 

WRITE  OUTPUT  TAPE  6.25  0178 

RETURN  0179 

END  0180 
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SUBROUTINE  GA I N ( KOP » DELTA » BETA »  FMC  »  S I GMA  »ER  »PHI »EL»H»RAIN)  0001 

C  POWER  GAIN  OF  ANTENNAS 

9749  FORMAT  (1F12.6)  0002 

RHI=PHI*. 01745329252  0003 

RELTA=DELTA*. 01745329252  0004 

GO  TO  (6*7*6) *KOP  0005 

6  T=COSF(RELTA)  0006 

RETA=BETA*. 01745329252  0007 

Q*SINF (RELTA)  0008 

R=Q*Q  0009 

S=R*R  0010 

X=18000**S IGMA/FMC  0011 

RHO=SQRTF( (ER-T*T)*(ER-T*T)+X*X)  0012 

RHOl 2=SQRTF ( RHO )  0013 

ALPHA=-ATANF (X/ ( ER-T*T ) )  0014 

PSIH=ATANF(2.*RH012*0*SINF( ALPHA*. 5 ) / ( RHO-R ) ;  0015 

CH=SQRT  F ( RHO*RHO+S-2 .*RHO*R*COSF ( ALPHA ) ) / ( RHO+R+2 .*RH012*Q*C0SF ( AL  0016 

1PHA**5>)  0017 

WAVE=299.7925/FMC  0018 

FAC =(3.1415926 54 *EL)/WAVE  0019 

U1=1.-T*SINF (RHI+RETA)  0020 

U2=1.-T*SINF (RHI-RETA)  0021 

EFF=0.  0022 

X=1 .  0023 

GO  TO  (2.2*3) .KOP  0024 

2  RAIN=( 3.2*COSF(RHI ) *COSF ( RHI ) *S I NF ( F AC*U1 ) *S I NF ( FAC*U 1 ) *S I NF ( FAC*U  0025 

12)*SINF(FAC*U2)*(CH»CH+l.-2.*CH*COSF(PS!H-(  12. 56637062‘*H)  /WAVE*Q)  )  0026 

2  )/(Ul*Ul*U2*U2)*(COSF(RETA)-SINF(RHI ) *T ) * ( COSF ( RET A ) “SI  NF ( RH I) *T )  0027 

ROK  =  3 »  0028 

GO  TO  4  0029 

3  RAI N= ( .5562474*COSF(1.570796*T)*COSF( 1.570796*T)*(CH*CH+1.-2.*CH*C  0030 

10SF ( PSIH-( 12  »  56637062 *H ) /WAVE*Q ) ) )/S  0031 

R0K  =  -4*  77  0032 

GO  TO  4  0033 

7  E=ER  0034 

I F ( H )  24.24*20  0035 

20  WAVE =29 9 *792 5 /FMC  0036 

A= (6.283185*H) /WAVE  0037 

X=H/WAVE  0038 

GO  TO  25  0039 
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24  EN  =  I  NTF (ABSF ( H ) / 1 0* )  0040 

IF(EN)  50 »  5 1  1 50  0041 

51  H1=ABSF ( H )  0042 

GO  TO  52  0043 

50  H1=EN/ ( ABSF ( H )  -lO.*EN)  0044 

52  A=H1*6. 283185  0045 

X  =H 1  0046 

25  D=2.*A  0047 

ROK=+5 •  0048 

2=2. *D  0049 

W=CI(Z»W1)  0050 

W2  =  C I ( D  »W3  )  0051 

R'V=30.*(-.5*COSF  (D)*(  »5772156649  +  LOGF(2)-W  )  +  (  l.+COSF(D)  >*<  .577  0052 

12156649+LOGF (D)-W2  ) +S I NF ( D ) * ( • 5*W 1 -W3 ) )  0053 

26  S=SINF (RELTA)  0054 

C=COSF ( REL  TA )  0055 

HM=A/6. 283185307  0056 

ANUM=COSF(A*S)-COSF(A)  0057 

28  CALL  VREFCO(RELTA»E.FMC»SIGMA.CV»PSlV)  0058 

FAC1=CV*CV  +  1 .-2.*CV*COSF ( PS  I V- 1 2 . 566  3 7  06*HM*S )  0059 

FAC2=2.*COSF( 6.2831853072*HM*S )  0060 

29  RAIN  =( 120.*ANUM*ANUM)/(RA*C*C)  0061 

RAIN=RAIN*FAC1  0062 

4  IF(RAIN)  11.11.31  0063 

31  R A I N  =  4 . 3429448 19*LOGF(RAIN)  0064 

I F ( X- • 2  5 )  69,70*70  0065 

69  EFF=( ( ( 6416. 702573*X-6091. 333295 )*X+2 179. 890548 >*X-364; 81 73803 )*X+  0066 

125.64620146  0067 

GO  TO  71  0068 

70  EFF=0.  0069 

71  RAIN=RAI N-EFF-ROK  0070 

41  IFIRAIN+10.  )11,11,10  0071 

11  RAIN =- 10.  0072 

10  RETURN  0073 

END  0074 
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SUBROUTINE  VREFCO ( DEL T A  * ER »FREQ »S I GMA »CV » PS  I V )  0001 

C  VERTICAL  GROUND  REFLECTION  COEFFICIENT 

PI=3. 141592654  0002 

X= 18 000. *S I GMA /FREQ  0003 

U=(ER*ER+X*X )  0004 

V=SQRTF ( U  )  0005 

Q=S I NF ( DEL  T  A )  0006 

R=Q*Q  0007 

S*R*R  0008 

T*COSF (DELTA)  0009 

RHO=SQR  TF ( ( ER-T*T ) * ( ER-T*T ) +X*X )  00 10 

RH012=SQRTF( RHO)  0011 

ALPHA=~ATANF(X/< ER-T*T ) )  0012 

A  =  2  • *RHOl2*Q*V*S I NF ( ALPHA*. 5+ AS  I NF ( X/V ) )  0013 

B=RHO-U*R  0014 

I F ( B )  10,20.30  0015 

20  I F ( A )  21,22,23  0016 

21  PSIV=-.5*PI  0017 

GO  TO  50  0018 

22  PSIV=0.  0019 

GO  TO  50  0020 

23  PSI V= • 5*P I  0021 

GO  TO  50  0022 

10  PSI V=ATANF ( A/B ) +P I  0023 

GO  TO  50  0024 

30  PSIV=ATANF(A/8)  0025 

50  CV  =  SQRTF(RH0*RH0+U*U*S-2.*RH0*U*R*C0SF(ALPHA+2.*ASINF(>X/V)  »  )  /  (  RHO+  0026 
1U*R+2.*RH012*V*Q*C0SF( ALPHA*. 5+AS INF (X/V) ) )  0027 

RETURN  0028 

END  0029 
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FUNCTION  C  I  <  X * S I  )  0001 

C  SINE  AND  COSINE  INTEGRAL 

4  IFIX-10. )5»60,60  0002 

5  SQ=X*X  0003 

CI=.5772156649+L0GF(X)  0004 

TERM*-l.*SQ/4.  0005 

G=4.  0006 

10  C I  =  C I +  T ERM  0007 

TERM=-l.*TERM*SQ*(G-2. )/< (G-1.)*G*G)  0008 

G=G+2 •  0009 

IFIABSF(TERM)-. 00005)  20,20»10  0010 

20  SI=0.  0011 

TERM=X  0012 

G=3.  0013 

25  S I =S I +TERM  0014 

TERM=-l.*TERM*SQ*(G-2. ) /( ( G— 1 • )*G*G)  0015 

G=G+2 •  0016 

IFIABSF(TERM)-. 00005)  80,80,25  0017 

60  TERM1=1.  0018 

T«l.  0019 

C=1 .  0020 

61  FAC=(4.*T*T-2.*T)/(X*X)  0021 

TERM2=-1.*FAC*TERM1  0022 

IF( ABSF(TERM2)-ABSF( TERM1 ) ) 62 *64,64  0023 

62  C=C+TERM2  0024 

T=T+1.  0025 

I F  < ABSF(TERM2/C)-10.E-10)65.65.6  3  0026 

63  TERM1=TERM2  0027 

GO  TO  61  0028 

64  SC=2.*T  0029 

THET  A  =  X-SC  0030 

PI =.5- < 1.-2«*ThETA)/(4.*SC)+( 1 . -4 . *ThETA-2 . *ThE TA*THET A ) / ( 8 . *SC*SC  0031 

1 )+( 3.+18.*THETA+8.*THETA*THETA) / ( 16.*SC*SC*SC )  0032 

C=C+TERM2*PI  0033 

65  TERM1=1./X  0034 

D=TERM1  0035 

T  =  1 .  0036 

66  FAC=(4,*T*T+2.*T)/(X*X)  0037 

TERM2=-1.*FAC*TERM1  0038 

IF(ABSF(TERM2)-ABSF(TERM1 ) >67,69,69  0039 
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67  D=D+TERM2  0040 

T=T+1 •  0041 

IF(ABSF(TERM2/D)-10.E-10) 70 *70 *68  0042 

68  TERM1=TERM2  0043 

GO  TO  66  0044 

69  SD=2.*T+1.  0045 

THET  A  =  X-SD  0046 

P I  =  •  5- <  l.-2.*THETA)/(4.*SD)  +  a.-4.*THETA-2.*THETA*THETA>/(8.*SD*SO  _  0047 

l)+(3.+18.*THETA+8.*THETA*THETA)/( 16«*SD*SD*SD)  0048 

D=D+TERM2*PI  0049 

70  CX=COSF ( X )  0050 

SX=S I NF ( X )  0051 

CI=(C*SX-D*CX)/X  0052 

S I  =  1.5 70796327- ( C*CX+D*SX  )/X  0053 

80  RETURN  0054 

END  0055 
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SUBROUTINE  LOSS! RELTA t ER » FREQ .S I GMA . Y )  OOOl 

GROUND  REFLECTION  LOSS 

DELTA=RELTA*. 01745329  0002 

PI=3. 141592654  0003 

Xa18000»*S IGMA/FREQ  0004 

Us ( ER*ER+X*X )  0005 

V*SQRTF ( U  )  0006 

Q=SI NF ( DELTA )  0007 

R=G*Q  0008 

S=R*R  0009 

T  *COSF ( DELTA )  0010 

RHO=SQRTF< ( ER-T*T ) * ( ER-T*T ) +X*X  )  0011 

RH012=SQRTF ( RHO)  0012 

ALPHA*-ATANF(X/I ER-T*T ) )  0013 

CH=SQRTF(RH0*RH0+S-2.*RH0*R*C0SF( ALPHA) ) / ( RHO+R+2 • *RH01 2*Q*C0SF ( AL  0014 
1PHA*#5  )  )  0015 

50  CV=SQRTF(RH0*RH0+U*U*S-2.*RH0*U*R*C0SFI ALPHA+2 . *AS I NF ( X/V ) ) ) / ( RHO+  0016 
1U*R+2.*RH012*V*Q*C0SF< ALPHA*. 5+ASINF(X/V) ) )  0017 

Y310.*LOG10F ( ( CH*CH+CV»CV ) /2. )  0018 

RETURN  0019 

END  0020 
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SUBROUTINE  POLYntJ,NN,MM*V»X»Y ) 

POWER  SERIES  VARIABLES 

0001 

DIMENSION  AU0*7,14) 

0002 

DIMENSION  P ( 29  » 16  *6 ) * ABP (2*6) 

0003 

COMMON  A8P.P 

0004 

COMMON  A 

0005 

Y*0 

0006 

M*MM 

0007 

120 

N=NN 

0008 

COEF*0 

0009 

127 

COEF«COEF*V+A(M»N*KJ) 

0010 

N=N-1 

oon 

I F ( N )  122*122,127 

0012 

122 

Y=Y*X+COEF 

0013 

M*M— 1 

0014 

I F ( M )  126*126,120 

0015 

126 

RETURN 

0016 

END 

0017 

-  89  - 


SUBROUTINE  NO  I SY  (  K.J  ,  XP  ,  CEG  .  ATNO  ) 

FOURIER  VARIABLES  AND  ATMOSPHERIC  RADIO  NOISE 

0001 

DIMENSION  P( 29*16*6) *ABP<  2*6 ) . ZZ ( 2 9 ) 

0002 

COMMON  ABP*P 

0003 

ALF=ABP ( 1 *KJ ) 

0004 

BET  =  ABP  <  2 . K J ) 

0005 

0=0* 008726646646 26*CEG 

0006 

C 1=C0SF { Q ) 

0007 

S 1=S I NF ( 0 ) 

0008 

DO  56  J  =  1 . 29 

0009 

• 

o 

II 

<* 

0010 

SX  =  S1 

0011 

CX  =  C1 

0012 

DO  55  K  =  1 *  15 

0013 

R=R+SX*P(J*K*KJ) 

0014 

SS=SX*C1+CX*S1 

0015 

CX=CX*C1-SX*S1 

0016 

sx=ss 

0017 

ZZ( J)=R+P( J.16»KJ) 

0018 

Q=. 01745329252*1 XP+90. ) 

0019 

S  1  =  S I NF ( 0 ) 

0020 

SX  =  S1 

0021 

Cl=COSF(Q) 

00  22 

CX=C1 

0023 

R  =  0. 

0024 

DO  57  K= 1 ♦ 29 

0025 

R  =R  +  SX*ZZ ( K ) 

0026 

SS=SX*C1+CX*S1 

0027 

CX=CX*C1-SX*S1 

0028 

SX  =  SS 

0029 

ATNO=R+ALF+BET*Q 

0030 

RETURN 

0031 

END 

0032 
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IX.  MATHEMATICAL  EXPRESSIONS 

1.  Great  circle  distance  in  kilometers,  degrees  and 
nautical  miles 

cos  (gcd)  =  sin  (xl)  •  sin  (x2)  +  cos  (xl) 
cos  (x2)  •  cos  (y2  -  yl) 
gcd  =  great  circle  distance  -  degrees 
xl  =  transmitter  latitude  -  degrees 
yl  =  transmitter  longitude  -  degrees 
x2  =  receiver  latitude  -  degrees 
y2  =  receiver  longitude  -  degrees 

2.  Bearing  from  receiver  to  transmitter,  degrees 
east  of  north.  (0  —  btry  ^  360°). 


btry  =  114.  5816  •  tan  * 


sin  (u-90  +  x2)  •  sin  (u-gcd) 


/  cind 


u  =  (180.  -  xl  -  x2  +  gcd)/  2 
cind  =  sin  (u)  •  sin  (u  -  90  +  xl) 

3.  Geographic  latitude  of  control  points  along  great 
circle  route. 

x'  1  =  cos  (pp)  •  cos  (90  -  x2)  +  sin  (pp)  •  sin  (90  -  x2) 

•  cos  (btry) 

(Control  Lat)  =  90  -  arccos  (x'  1) 

pp  =  great  circle  distance  in  degrees  from  transmitter 
terminal  to  control  point 

4.  Geographic  longitude  of  control  points  along  great 
circle  route,  (y  =  y2  -  y'  ) 
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y' 


arccos 


cos  (pp)  -  cos  (90  -  x2) 


•  cos  (90  -  x3)/ 

sin  (90  -  x2)  •  sin  (90  -  x3) 


x3  =  control  point  geographic  latitude 

5.  Geomagnetic  latitude  of  control  points  along  great 
circle  route,  (y  =  90  -  y'  ) 


y' 


arccos 


sin  (x4)  •  sin  (x3)  + 
cos  (x4)  •  cos  (x3 ) 
cos  (x3  -  x4 ) 


x4  =  latitude  of  geomagnetic  north  pole 

6.  Local  time  at  receiver  terminal. 

LMT  =  GMT  -  y2/l  5 

LMT  =  local  mean  time  -  hours 
GMT  -  Greenwich  mean  time  -  hours 

7.  Sun' s  zenith  angle  at  control  points. 

cos  (ijj)  =  sin  (z)  .  sin  (ssp)  +  cos  (z)  .  cos  (ssp)  . 
cos  [  (15  x  GMT)  -  180]. 

4<  =  sun' s  zenith  angle  -  degrees 

z  =  control  point  latitude  -  degrees 

ssp  =  latitude  of  subsolar  point  of  sun  for  middle 

of  month  in  question  -  degrees 

GMT  -  Greenwich  mean  time  -  hours 

8.  Ionospheric  absorption  index  "I" 

I  =  (1  +  .  0037  SSN)  (cos  .  881  3 


-  92  - 


9. 


SSN  4»  predicted  or  observed  12  month  moving  average 
Zurich  sunspot  number.  (0  ~  SSN  ^  200) 


E -layer  distance  factor. 


2. 085000000  • 

10 

~l  -  1.  33200756  •  lO-2  * 

x  + 

-22  ... 

-3 

3 

3. 761385053  • 

10 

x  -  5.  038476266  •  10 

x  + 

-4  4 

-6 

5 

2. 624808315  * 

10 

x  -  5.  976618436  *  10 

x  + 

1. 334494261  • 

10 

7  6 

x  -  4.  368460907  *  10 

-9 

• 

7 

X 

x  =  great  circle  distance  -  hundreds  of  statute  miles 

10.  E -2000  MUF. 

emc  =  3.  345996232  +  37.  67736072  •  1  -  52.41191754*  I2 
+  39.26151056  •  I3  -  10.66484988  *  I4 


I  =  absorption  index 

11.  E-layer  MUF. 
emufy  =  (emc)  •  (elfc) 
emc  =  E-2000  MUF 

12.  F-layer  distance  factor. 

fife  =  4.  699243101  •  10~3  x  +  2.  264634341  *  10"3  *  x2 

+  9.202437332  *  10“5  *  x3  +  6.  86525981 7  *  10 “5  *  x4  - 
9.985831104  •  10"6  •  x5  +  4.  49 1 51441  10”7  •  x6  - 

6.  712654756  •  10“9  *  x? 

x  =  great  circle  distance  -  hundreds  of  statute  miles 

13.  F 2 -layer  Fourier  generation  of  foF2  and  M-3000  factors. 
Fourier  time  variation  function,  used  to  obtain  foF2 

or  M-3000  factor: 
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J 


(x»  y 


jb 


(x,y)  cos  (jb  -  1)  t  + 


jb  =  2 


ba.,  (x,y)  sin  (jb  -  1)  t 
Jb 


I 


is  =  2jb  -  1 
jb  =  1, 2, 3.  ... j 


ka  =  1 


I 


Q 


is ,  ka,  io  ka 


(xy)  is  =  2jb  -  2 


ka  =  1 


jb  -  2 ,  3 , 4.  .  .  .  j 


x  =  geographic  latitude  of  control  point  -  degrees 
(-90°  <  x  ^  90°) 

y  =  geographic  longitude  of  control  point  -  degrees 
(0°  <  x  ^  360°  -  East  of  North) 

t  =  local  time  at  receiver  -  (hour  angle)  (180  —  t  —  180°) 
I,j,k,L  =  constants  describing  number  of  harmonics 
in  the  Fourier  functions 

io  =  matrix  for  foF2  or  M-3000  (high  or  low  SSN) 

Q.  =  Fourier  coefficients  defining  the  function 

is,  ka, io 


(x,y,t)  (Time  variation) 


Latitudinal  and  longitudinal  variation  functions  are 
shown  in  Tables  1  and  2. 


For  a  detailed  explanation  of  the  generation  of  (x,y,t) 


refer  to  [Jones  and  Gallet,  1961].  The  above  formulas  and 
constants  have  been  altered  from  the  originals  in  the  reference 
to  make  the  generation  permissible  on  a  larger  variety  of 
computers . 
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Table  1.  Geographical  Function  in  Latitude 


MAIN  LATITUDINAL  VARIATION 

ka 

Gka<x' y> 

1 

1 

2 

sin  (x) 

3 

sin^(x) 

*  *  * 

k 

sin  (x) 
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Table  2.  Geographical  Functions  in  Latitude 
and  Longitude. 


MIXED  LATITUDINAL  AND  LONGITUDINAL 
VARIATION  -  FIRST  ORDER  IN  LONGITUDE 


ka 

Gka(x’y) 

k  +  1 

cos  (x)  cos  (y) 

k  +  2 

cos  (x)  sin  (y) 

k  +  3 

sin  (x)  cos  (x)  cos  (y) 

k  +  4 

sin  (x)  cos  (x)  sin  (y) 

.  .  . 

L  -  1 

sin  (x)  cos  (x)  cos  (y) 

L 

sin  (x)  cos  (x)  sin  (y) 

X 

i 

L 

n 

nj 

- ^ -  -  1 

2 

MIXED  LATITUDINAL  AND  LONGITUDINAL 
VARIATION  -  SECOND  ORDER  IN  LONGITUDE 


2 

L  +  1 

cos  (x)  cos  (2  y) 

L  +  2 

2 

cos  (x)  sin  (2  y) 

L  +  3 

2 

sin  (x)  cos  (x)  cos  (2  y) 

L  +  4 

sin  (x)  cos^  (x)  sin  (2  y) 

.  .  . 

Lb  2 

I  -  1 

sm  (x)  cos  (x)  cos  (2  y) 

Lb  2 

I 

sm  (x)  cos  (x)  sin  (2  y) 

Lb  =  I  -  L 

— = -  -  1 

2 
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14. 


F2-layer  gyro  frequency  may  be  adequately  represented 
for  MUF -FOT  predictions  by  a  set  of  least  squares 
coefficients  describing  the  orthogonal  polynomial. 


y  =  (ai,i  + 


al,2  X 


ai,3  X 


ai ,  7  X  )  Z  +  •“ 


a7,  2  X 


+  a 


7,3 


2 

z 


7,7 


6 

x 


15. 

16. 


17. 


y  =  gyro  frequency  of  F2  layer  -  Me/  s 
x  =  longitude  of  control  point  -  degrees 
z  =  latitude  of  control  point  =  degrees 
a  =  set  of  least  squares  coefficients  describing  the 
orthogonal  polynomial  (y) 

F2-4000  MUF. 

F2-4000  MUF  =  (foF2)  x  (M-4000  factor) 


F -MUF  for  low  and  high  solar  activity. 

F -MUF  =  ZDF  +  fife  •  (F2-4000)  -  ZDF 

ZDF  =  foF2  +  \  F2-layer  gyro  frequency 
M-4000  factor  =  (1.  1)  *  M-3000  factor 
Interpolation  for  intermediate  values  of  solar  activity. 
(10  <  SSN  <  180) 


F -MUF  = 


F  -  MUF  •  (  180  -  SSN)  +  F-MUF^(SSN-10) 


/170 


F  -  MUF.  =  F -MUF  for  SSN  10 

(i) 

(H) 


F  -  MUF 


=  F -MUF  for  SSN  180 
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18.  Angle  at  the  ionosphere  (4>). 


4>'  =  sin 


-1 


- 1  A1 

p  (  sin  (90  +  tan  (cos  /  2  -  p/  (p  +  H  ) 


sin  A' 


12 


p  +  H 


S. 


4> '  =  angle  of  incidence  -  degrees 
p  =  radius  of  earth  -  km 
A1  =  great  circle  distance  -  degrees 
--  layer  height  -  km 


19. 


20. 


Ionospheric  absorption 
A  =  615.  5  N  sec  cj>  (I) 


(£  +  f  ) 
g 


1.  98 


cj>  =  angle  of  incidence  for  100  km  -  degrees 
f  =  operating  frequency 
I  =  absorption  index 
f  =  gyro  frequency  -  Me/  s 


Basic  transmission  loss  for  isotropic  antennas  in  free  space. 

Lbf  =  10  log10  (Pr/pJ  =  10  log10  (4  "  dA)2 


=  36.  58  +  20  log10(d)+20  l°g10(fMc/s> 


Pr  =  power  available  at  receiving  antenna 

p^  =  power  delivered  to  transmitting  antenna 

d  =  ray  path  distance  -  miles 
X.  =  wave  length 

fMc/  s  =  frecluency  "  Me  /  s 
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21.  Relationship  between  cj>'  and  A. 

sin  cj>'  =  p  cos  A 
P+H, 


4>'  =  angle  of  incidence  at  ionosphere 
p  =  radius  of  earth 
=  layer  height 

A  =  radiation  angle  of  wave 

22.  Ground  reflection  factors  for  vertical  and  horizontal 
polarization. 


sin  A 


sin  A  +  \[ 


k 


ix)  -  cos' 


-  ix)  -  cos 


A 


A 


(t  -  ix)  sin  A  -  \T,  T~  2 

..  r  (e  -  ix)  =  cos  A 

Ky  “  _ E _ 


(«  -  ix)  sin  A  +  \T  .  .  2 

r  (e  -  ix)  -  cos  A 

r 


€  =  relative  dielectric  constant  of  earth, 
r 

x  -  — - —  ~  18  x  103  (J/f, 

CO  €  — 

V 

O’  =  conductivity  of  earth  (mhos  /  meter), 
co  =  angular  frequency, 
f  =  frequency  in  megacycles, 

A  =  angle  of  elevation  in  degrees,  and 

i  =  \T  - 1. 

=  dielectric  constant  of  free  space 
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23.  Rhombic  antenna  power  gain  relative  to  isotropic  in  free  space. 


g(A,  P)  =  3.2  (HL-) 


{Dy(A,p)}2  +  {DH(A,p}2 


r-x  /a  o\  ,  sin  u  sin  u_  .  . 

Dy(A,  p)  =  cos  9  •  _ 1_  • _ 2_  •  sm  p  •  sin  A 

u  u 

1  2 


Kv|2  +1-2  |  Ky  |  cos  (4jv  -  •  sin  A) 


sin  sm  u^ 

D  (A, (3)  =  cos  4>  - - -  -  •  (cos  p  -  sin<j>  •  cos  A) 

H  U1  u2 


Kh|2  9  1  -  2|Kh|  cos  (^  -  ^  •  sin  A) 


Tri 

U1  "  X 


1  -  cos  A  •  sin  ( 4>  +  p) 


u 


Tri 

2  J  T 


1  -  cos  A  •  sin  -  p) 


g  (A,  p)  =  power  gain  relative  to  isotropic  in  free  space 
A  =  angle  of  departure  in  degrees, 

P  =  angle  of  azimuth  in  degrees, 
i  =  rhombic  leg  lenth  in  meters, 
h  =  antenna  height  in  meters, 

4>  =  tilt  angle  in  degrees. 


100  - 


r  2p; 


‘J'y  =  tan 


sin  A  .  sin  (tan 


-1  x 


+  -r> 


p  -  y  .  sin 


+  TT 


^  4  Z 

Y  |  _  [p  +  sin  A  -  2p  •  sin  A  •  cos  a  ] 

H1  =  r  .  2  A  71  ~  7~n 

[p  +  sm  A  +  ip  •  sm  A  •  cos  —  J 


K. 


2  2  .4  .2  .  —  .  -lx) 

p  +  y  sm  A  -  2p  .  y  .  sm  A  .  cos  (a  +  2  sm  — i 

_  y2 


.  2  A  _  2  T  .  A  .  c  .-lx. 

p+y.sin  A+2  y  .  sm  A  .  cos  (r  +  sm  -y-) 

P  2  2 

y 


=  tan 


-1 


2p  2  .  sin  A 


p  -  sm 


2  2 
y  =  e  +  x 
r 


p  =  [  {  e  -  cos2  A)2  +  x2] 
r 


1 

2 


-1 

a  =  -tan  x 

2 

(e  -  cos  A) 
r 

x  =  18  x  103  07  f 

€  =  relative  dielectric  constant  of  the  ground 

r  ° 

0"  =  conductivity  of  the  ground  in  mhos  per  meter 
f  =  operating  frequency  in  megacycles  per  second 

i  |  i  ili 

=  -  |  |  e  H  =  horizontal  reflection  coefficient,  and 
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Ky  =  -  |  Ky  |  e^V  =  vertical  reflection  coefficient. 


24.  Power  gain  of  half-wave  horizontal  dipole. 


=  (.  74582) 


cos  ( •  cos  A)  +  1-2  •  •  cos 


.  2  A 
sm  A 


(+■ 


4ttH 


H  \ 


sin  A) 


1 


K  =  amplitude  of  horizontal  reflection  coefficient 
ri 

ikj.  =  phase  amplitude  of  horizontal  reflection  coefficient 
ri 

H  =  height  of  antenna  -  meters 
\  =  299.  7925/frequency  in  megacycles 

H-* 


A  =  angle  of  elevation 

Note:  72  ohms  assumed  impedance  of  antenna 

Gain  =  10.  x  log10  (E  ) 

decibels  B10  A7 


2  5.  Power  gain  of  vertical  antennas 
120 


ea  ra 


cos  (A)  sin  A)  -  cos  (A) 
cos  (A) 


f  2 

lkv  + 


-  2  •  Ky  cos  (4<v  -  2A  •  sin  A  | 


V 


R  =30 
a 


cos  (2A) 


'  |c  +  In  (4A)  -  C.  (4A)| 


+  -I  I  +  cos  (2A)  1  x  'c  +  In  (2A)  -  C.  (2A) 

j  L  1  j 
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+  sin  (2A)  x 


Si  (4A)  -  S  (2A)  } 

2  J 


C  =  .  5772156649 
H  =  height  of  antenna  -  meters 
Si  =  sine  integral 
Ci  =  cosine  integral 
A  =  2  tt  H/\ 

26.  Efficiency  factor  for  short  vertical  antennas. 
EFF  =  25.  64620146  -  364.  8173803  x 

+  2179. 890548  x2  -  6091. 333295  x3 
+  6416. 702573  x4 
x  =  wave  length  of  antenna 
(1/16  ^  x  ==  1/4) 

27.  Ground  reflection  loss. 


G.  =  10 

f 


log 


10 


2  2 

kh  +kv 


28.  Relationship  of  field  strength  to  transmission  loss. 

(Decibels  above  isotropic) 

E  =  107.  2  +  20  log  .  _  (f  )  -  L. 

B10  '  me  b 

f  =  frequency  -  Me/  s 
me 

L,  =  transmission  loss  -  decibels 
b 
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X.  GENERALIZED  BLOCK  DIAGRAM  OF  HF  SYSTEM  PERFORMANCE  ROUTINE 


START 


Calculate 

1.  Great  circle  distance 

2.  Bearings 

3.  Location  of  control  points 

4.  LMT  at  control  points 

5.  Sun' s  zenith  angle  at  control 

points 

6.  Geomagnetic  location  of 

control  points 

7.  Absorption  index  I  at  control 

points 

8.  Path  absorption  index  I 

9.  E-layer  control  points  MUFs 

10.  F2-layer  control  points  MUFs 

11.  MUF  and  FOT  for  path 


C  alculate 

1.  F 2 -layer  height 

2.  Distance  per  hop 

3.  Radiation  angle  for  E,  EF2 

and  F2-layer  modes 

4.  Angle  of  incidence  associated 

with  mode 

5.  Maximum  E -layer  angle  of 

incidence 

6.  Minimum  F2- layer  angle  of 

incidence 

7.  Dominant  type  of  terrain 

8.  Possible  modes 

9.  Ionospheric  loss  for  path 

10.  Critical  frequency  E-layer 

11.  Critical  frequency  F2-layer 

12.  Ground  reflection  loss 

13.  Ray-path  attenuation 

14.  Transmitter  antenna  gain 

15.  Receiver  antenna  response 

16.  Mode  with  least  loss 


Read 

Graph  or  print 

1.  Circuit  input  parameters 

MUF -FOT  option 

2.  Ionospheric,  geophysical  and 

Figure  1  or  9 

systems  parameters 

"END" 

Print 

System  loss  option 
Figure  5 
"END" 


1.  Set  receiving  antenna 

response  to  isotropic 

2.  Relate  system  loss  to 

field  strength  (dbu) 

Print 

Field  strength 
option  Figure  6 
"END" 

Calcu  late 

1.  LMT  at  receiver  site 

2.  Atmospheric,  man-made 

and  galactic  noise 

3.  Compare  signal  and  noise 

Print 

S/N  ratio  option 
Figure  7 
"END" 

Calculate 

1.  Type  of  path,  e.  g.  ,  polar, 

short  or  long,  etc. 

2.  Circuit  reliability 

Print 

Reliability  option 

F igure  2  or  8 
"END" 

Calculate 

1.  Lowest  frequency  in  Band 

7  with  reliability  of  90% 

or  greater 

2.  Iterate  downward  on  freq. 

until  reliability  is  90% 

Graph  or  print 
FOT-LUF  option 
Figure  3,  4  or  10 
"END" 

Figure  18 
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XI.  CONCLUSIONS 

The  computer  solution  of  HF  systems  problems  based  on  estab¬ 
lished  manual  methods  is  efficient  and  practical. 

The  radio  systems  engineer  may  with  this  basic  tool  check  pre¬ 
dictions  against  operational  data  to  update  the  prediction  scheme  and 
to  test  new  prediction  parameters  with  little  effort. 

The  routine  is  based  on  average  monthly  values;  therefore,  it  is 
most  useful  in  systems  problems  such  as  allocation  of  frequencies  and 
circuit  design.  Use  of  the  relative  values  produced  by  the  prediction 
scheme  are  no  doubt  more  important  to  the  communication  engineer 
than  are  absolute  values.  The  routine  is  most  valuable  for  the  exper¬ 
ienced  engineer  with  an  adequate  knowledge  of  the  shortcomings,  as 
well  as  the  usefulness  of  such  a  prediction  scheme. 
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In  general,  each  section  carries  out  specialized  research,  development,  and  engineering  in  the  field  indicated  by 
its  title.  A  brief  description  of  the  activities,  and  of  the  resultant  publications,  appears  on  the  inside  of  the 
front  cover. 


WASHINGTON,  D.C. 

Electricity.  Resistance  and  Reactance.  Electrochemistry.  Electrical  Instruments.  Magnetic  Measurements. 
Dielectrics.  High  Voltage.  Absolute  Electrical  Measurements. 

Metrology.  Photometry  and  Colorimetry.  Refractometry.  Photographic  Research.  Length.  Engineering  Metrology. 
Mass  and  Volume. 

Heat.  Temperature  Physics.  Heat  Measurements.  Cryogenic  Physics.  Equation  of  State.  Statistical  Physics. 
Radiation  Physics.  X-ray.  Radioactivity.  Radiation  Theory.  High  Energy  Radiation.  Radiological  Equipment. 
Nucleonic  Instrumentation.  Neutron  Physics. 

Analytical  and  Inorganic  Chemistry.  Pure  Substances.  Spectrochemistry.  Solution  Chemistry.  Standard  Refer¬ 
ence  Materials.  Applied  Analytical  Research.  Crystal  Chemistry. 

Mechanics.  Sound.  Pressure  and  Vacuum.  Fluid  Mechanics..  Engineering  Mechanics.  Rheology.  Combustion 
Controls. 

Polymers.  Macromolecules:  Synthesis  and  Structure.  Polymer  Chemistry.  Polymer  Physics.  Polymer  Charac¬ 
terization.  Polymer  Evaluation  and  Testing.  Applied  Polymer  Standards  and  Research.  Dental  Research. 
Metallurgy.  Engineering  Metallurgy.  Metal  Reactions.  Metal  Physics.  Electrolysis  and  Metal  Deposition. 
Inorganic  Solids.  Engineering  Ceramics.  Glass.  Solid  State  Chemistry.  Crystal  Growth.  Physical  Properties. 
Crystallography. 

Building  Research.  Structural  Engineering.  Fire  Research.  Mechanical  Systems.  Organic  Building  Materials. 
Codes  and  Safety  Standards.  Heat  Transfer.  Inorganic  Building  Materials.  Metallic  Building  Materials. 

Applied  Mathematics.  Numerical  Analysis.  Computation.  Statistical  Engineering.  Mathematical  Physics.  Op¬ 
erations  Research. 

Data  Processing  Systems.  Components  and  Techniques.  Computer  Technology.  Measurements  Automation. 
Engineering  Applications.  Systems  Analysis. 

Atomic  Physics.  Spectroscopy.  Infrared  Spectroscopy.  Far  Ultraviolet  Physics.  Solid  State  Physics.  Electron 
Physics.  Atomic  Physics.  Plasma  Spectroscopy. 

Instrumentation.  Engineering  Electronics.  Electron  Devices.  Electronic  Instrumentation.  Mechanical  Instru¬ 
ments.  Basic  Instrumentation. 

Physical  Chemistry.  Thermochemistry.  Surface  Chemistry.  Organic  Chemistry.  Molecular  Spectroscopy.  Ele¬ 
mentary  Processes.  Mass  Spectrometry.  Photochemistry  and  Radiation  Chemistry. 

Office  of  Weights  and  Measures. 

BOLLDEK,  COLO. 

CRYOGENIC  ENGINEERING  LABORATORY 

Cryogenic*  Processes.  Cryogenic  Properties  of  Solids.  Cryogenic  Technical  Services.  Properties  of  Cryogenic 
Fluids. 

CENTRAL  RADIO  PROPAGATION  LABORATORY 

Ionosphere  Research  and  Propagation.  Low  Frequency  and  Very  Low  Frequency  Research.  Ionosphere  Re¬ 
search.  Prediction  Services.  Sun-Earth  Relationships.  Field  Engineering.  Radio  Warning  Services.  Vertical 
Soundings  Research. 

Troposphere  and  Space  Telecommunications.  Data  Reduction  Instrumentation.  Radio  Noise.  Tropospheric 
Measurements.  Tropospheric  Analysis.  Spectrum  Utilization  Research.  Radio-Meteorology.  Lower  Atmosphere 
Physics. 

Radio  Systems.  Applied  Electromagnetic  Theory.  High  Frequency  and  Very  High  Frequency  Research.  Fre¬ 
quency  Utilization.  Modulation  Research.  Antenna  Research.  Radiodetermination. 

Upper  Atmosphere  and  Space  Physics.  Upper  Atmosphere  and  Plasma  Physics.  High  Latitude  Ionosphere 
Physics.  Ionosphere  and  Exosphere  Scatter.  Airglow  and  Aurora.  Ionospheric  Radio  Astronomy. 

RADIO  STANDARDS  LABORATORY 

Radio  Standards  Physics.  Frequency  and  Time  Disseminations.  Radio  and  Microwave  Materials.  Atomic  Fre¬ 
quency  and  Time-Interval  Standards.  Radio  Plasma.  Microwave  Physics. 

Radio  Standards  Engineering.  High  Frequency  Electrical  Standards.  High  Frequency  Calibration  Services.  High 
Frequency  Impedance  Standards.  Microwave  Calibration  Services.  Microwave  Circuit  Standards.  Low  Frequency 
Calibration  Services. 

Joint  Institute  for  Laboratory  Astrophysics- NBS  Group  (Univ.  of  Colo.). 


